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Abstract 

The  one  dimensional  initial  value  problem  of  a  monatomic 
single  component  gas  is  considered.   Using  the  linearized  Boltzmann 
equation  the  dispersion  relation  is  studied.   In  addition  to  the 
usual  gas  dynamic  sound  waves  one  finds  an  infinity  of  decaying 
propagating  waves.   The  phenomenon  naturally  exhibits  itself  as 
a  sequence  of  epochs,  the  last  stage  of  which  is  hydrodynamic . 

With  reference  to  the  same  problem  macroscopic  equations 
such  as  Euler,  Navier-Stokes,  Burnett,  Grad '  s  moments  equations, 
etc.,  are  considered.   In  addition  the  recently  considered 
"kinetic  models"  of  Gross  et  al .  are  applied  to  the  problem. 
These  various  formulations  are  critically  analyzed  and  compared 
with  each  other  and  with  the  Boltzmann  analysis.   Lastly,  several 
alternate  molecular  and  macroscopic  equations  are  offered  which 
remedy  some  of  the  shortcomings  which  appear  in  the  above 
mentioned  approximate  theories. 
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The  Initial  Value  Problem,  Sound  Propagation,  and 
Modeling  in  Kinetic  Theory 

1 .   Introduction 

The  mathematical  difficulties  encountered  in  dealing  with 
the  Boltzmann  equation  have  motivated  the  creation  and  study  of 
many  approximate  formulations.   In  this  class  of  approximate 

formulations  are  the  Euler  equations,  Navier-Stokes  equations, 

1  2 

Burnett  equations  ,  the  moments  equations  of  H.  Grad  ,  and  the 

interpolated  equations  of  H.  Grad.    A  feature  shared  by  each 

of  these  systems  is  that  they  may  be  considered  as  having  been 

extracted  from  the  Boltzmann  equation  by  a  process  which  regards 

a  flow  as  being  smooth  with  respect  to  the  mean-free -path  as  a 

length  scale,  and  smooth  with  respect  to  the  collision  frequency 


1.  S.  Chapman  and  T.  G.  Cowling,  The  Mathematical  Theory  of 
Non -Uniform  Gases,  Cambridge  University  Press  (1952) . 

A  fluid  dynamicist  might  object  to  classifying  the  Navier- 
Stokes  equations  as  an  approximation  to  the  Boltzmann  equation. 
His  objection  is  warranted  by  the  fact  that  the  Navier-Stokes 
equations  apply  to  a  wider  range  of  fluids  than  does  the 
Boltzmann  equation.   We,  however,  restrict  our  attantion  to 
gases,  in  which  case  the  Euler  equations,  Navier-Stokes  equa- 
tions, and  Burnett  equations  can  be  thought  of  as  being 
extracted  from  the  Boltzmann  equation  by  the  Chapman -Enskog 
procedure  as  given,  for  instance,  in  the  cited  reference. 

2.  H.  Grad,  Comm.  Pure  and  Appl .  Math.  2,  331  (19^9)- 

3-   H.  Grad,  "Theory  of  Rarefied  Gases",  in  Rarefied  Gas  Dynamics, 
ed.  F.  M.  Devienne  (Pergamon  Press,  London,  I960) ,  p.  100. 

Grad  gives  a  technique  whereby  the  thirteen  moments  equations 
in  a  quick  way  give  the  Navier-Stokes  equation.   This  method, 
which  Grad  calls  interpolation,  is  generalized  by  him  to  any 
system  of  monents  equations. 
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as  a  time  scale.   This  presumably  excludes  from  consideration 
such  "kinetic"  boundary  layers  as  the  region  adjacent  to  a  wall 
(not  to  be  confused  with  the  fluid  mechanical  boundary  layer 
which,  in  general,  is  much  thicker),  the  interior  of  a  strong 
shock,  and  layers  of  large  shear.   There  is,  in  addition  to  the 
spatial  boundary  layers  just  mentioned,  the  associated  temporal 
boundary  layer.   This  boundary  layer  arises  from  the  fact  that 
the  above  mentioned  techniques  for  deriving  fluid  equations  are 
presumably  not  applicable  at  all  for  times  less  than  the  mean 
collision  time. 

An  entirely  different  class  of  approximate  formulations 
has  its  origin  in  the  model  equation  associated  with  Bhatnager, 

Gross,  and  Krook  ,  and  Welander.    This  was  subsequently  gener- 

6 
alized  by  Gross  and  Jackson  to  give  a  wide  range  of  model 

equations.   Whereas  the  approximate  formulations  are  equations 
in  moments  of  the  molecular  distribution  function,  a  model  equa- 
tion is  an  equation  in  the  molecular  distribution  function 
itself.   These  model  equations  share  many  common  features  with 
the  Boltzmann  equation,  and  are  sometimes  more  amenable  to  analysis 
As  will  be  shown,  these  model  equations  overcome  some  serious 
shortcomings  of  the  macroscopic  equations  (and  introduce  others). 


4.   P.  P.  Bhatnager,  E.  P.  Gross,  and  M.  Krook,  Phys.  Rev.  £4, 
511  (195^). 

5-   P.  Welander,  Arkiv  Fysik  J_,  No.  44  (1954). 

6.   E.  P.  Gross  and  E.  A.  Jackson,  Phys.  Fluids  2,  4^2  (1959)- 
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One  of  the  major  objectives  of  this  study  is  to  relate 
the  above  mentioned  approximations  (and  others)  to  the  Boltzmann 
equations  and  to  each  other.   A  precise  analysis  would  involve 
estimates  on  the  solutions  of  all  systems  under  consideration 
and  a  subsequent  comparison  of  them.   For  obvious  reasons  such  a 
plan  is  not  feasible.   Instead  we  will  attempt  a  more  modest  and 
in  some  respects  a  more  revealing  approach.   The  non-linear  theory 
will  be  discarded  and  we  consider  only  the  initial  value  problem 
of  the  one  dimensional  linearized  theory  in  the  absence  of 
boundaries.   Our  analysis  will  be  performed  exclusively  for  this 
elementary  but  by  no  means  simple  problem.   We  will  consider  the 
latter  problem  as  exhaustively  as  possible,  and  this  will  take 
up  a  considerable  part  of  this  report.   The  approximate  theories 
will  then  be  applied  to  the  same  initial  value  problem.   Only 
when  all  of  these  calculations  have  been  performed  will  we  be 
able  to  connect  the  various  theories  amongst  themselves  and  to  the 
Boltzmann  equation.   As  a  word  of  caution  it  should  be  noted  that 
the  various  interrelations  which  are  made  can  only  justifiably  be 
made  for  the  particular  case  under  study  --  it  is  tempting  none- 
theless to  consider  these  results  as  being  more  general.   On 
the  other  hand,  any  inadequacy  of  a  set  of  equations  which  is 
pointed  out  is  an  inadequacy  in  general,  although  it  may  be  less 
important  In  another  problem. 

The  analysis  of  the  Boltzmann  equation  has  been  based  on  the 
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investigations  of  sound  propagation  made  by  C.  S.  Wang  Chang  and 

7 
G.  E.  Uhlenbeck  ,  and  the  later,  more  general  developments  in 

o 

linearized  kinetic  theory  made  by  H.  N.  Mott -Smith.    The  author 
is  in  no  small  way  indebted  to  these  workers  as  is  evidenced  by 
the  frequent  reference  to  their  work. 

Before  going  on  to  the  analysis  it  will  be  advantageous  to 
outline  some  of  the  results  which  are  later  found.   The  evolution 
in  time  of  the  solution  from  the  initial  time  is  in  a  natural 
manner  divided  into  epochs.   With  each  epoch  is  associated  a 
propagating  mode  or  a  diffusing  mode  or  both,  there  being  an  infinity 
of  each  type.   As  time  evolves,  clusters  of  modes  become  washed 
out  of  the  picture  in  the  sense  that  they  have  a  shorter  exponen- 
tial damping  time  than  those  remaining .   In  the  final  stage  there 
is  just  one  cluster  of  modes  compared  to  which  all  other  modes 
are  asymptotically  small.   This  is  the  usual  hydrodynamic  regime. 

Macroscopic  equations  such  as  the  Navier-Stokes  equations 
and  Burnett  equations  only  furnish  a  description  of  the  hydro- 
dynamic  stage  of  the  initial  value  problem.   Higher  approximations 
such  as  Burnett  and  any  other  that  might  come  from  a  Chapman- 
Enskog  procedure  attempt  to  give  a  finer  description  of  this  last 

stage . 

2 
The  moments  equations  of  H.  Grad  ,  on  the  other  hand,  attempt 


7.  C.  S.  Wang  Chang  and  G.  E.  Uhlenbeck,  "On  the  Propagation  of 
Sound  in  a  Monatomic  Gas",  Univ.  of  Mich.,  Engr.  Rep.  (1952). 

In  this  highly  original  and  stimulating  report,  the  authors 
make  the  first  general  attack  on  sound  propagation,  by  means 
of  the  Boltzmann  equation.   The  analysis  in  Section  2  is 
largely  based  on  the  work  of  this  report,  and  that  of  the 
following  cited  work  of  H.  M.  Mot t -Smith. 

8.  H.  M.  Mott -Smith,  "A  New  Approach  to  the  Kinetic  Theory  of 
Gases",  Lincoln  Laboratory,  MIT  (195^)- 
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to  take  into  consideration  some  of  the  modes  that  do  occur  prior 
to  the  hydrodynamic  stage.   However  these  equations  have  two 
shortcomings .  -  First  of  all  any  moments  system  has  a  bound  on 
the  absorption  of  the  system.   This  is  contrary  to  the  results 
gotten  from  Boltzmann  equation,  which  has  modes  of  arbitrary  high 
absorption,  and  suggests  moreover  that  any  mode  approaches 
unbounded  absorption  with  high  frequencies.   As  is  shown  the 
interpolation  scheme  of  H.  Grad   overcomes  this  by  changing  the 
type  of  the  equation.   For  instance  interpolation  sends  the 
thirteen  moments  equations  into  the  Navier-Stokes  equations. 
The  second  drawback  to  the  moment  methods  involves  the  onset  of 
propagation  for  the  non-hydrodynamic  mode.   They  can  incorrectly 
predict  a  threshold  frequency  at  which  extraordinary  propagation 
occurs.   The  latter  effect  is  connected  to  the  moments  one 
chooses  in  a  description.   Alternative  systems  of  equations  can 
be  introduced  which  overcome  these  difficulties. 

The  equations  falling  under  the  title  of  kinetic  models 
also  take  into  consideration  the  presence  of  other  modes.   It 
does  this  by  lumping  together  an  infinity  of  modes  into  a  single 
cluster  and  giving  a  more  accurate  description  of  a  few  modes. 
The  single  relaxation  model  of  Bhatnager,  Gross  and  Krook  for 
instance  gives  only  a  qualitatively  correct  description  of  the 
hydrodynamic  stage.   The  kinetic  models  which  occur  in  the 
literature  give  an  incorrect  picture  of  the  onset  of  the 
extraordinary  sound  propagations.   They  do  this  for  essentially 
the  same  reason  as  this  occurs  in  the  moments  equations. 


Alternate  models  are  easily  developed  which  overcome  this,  and 
without  a  loss  of  simplicity.    Moreover  kinetic  models  also 
share  with  the  moments  systems  the  property  of  placing  a  bound 
on  the  absorption  of  the  system.   On  the  other  hand  within  the 
confines  of  this  absorption  band  the  kinetic  models  give  at 
least  a  qualitatively  more  correct  description  of  the  spectrum 
of  the  linearized  Boltzmann  equation.   Furthermore  kinetic 
models  can  be  applied  to  highly  rarefied  gases.   For  instance 
they  can  be  used  in  a  Knudsen  iteration  procedure.   Macroscopic 
equations  by  their  very  nature  cannot  be  so  applied. 

In  regards  to  the  initial  value  problem,  kinetic  models 
lie  between  the  Boltzmann  equation  analysis  and  a  macroscopic 
description.   In  fact  by  choosing  a  sufficiently  detailed 
kinetic  model  any  macroscopic  system  of  equations  can  be 
extracted  from  it  without  any  resort  to  the  Boltzmann  equation. 

Acknowledgment:   The  author  is  grateful  to  C.  A.  Sirovich  for 
many  helpful  mathematical  discussions. 
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2.  The  Linearized  Boltzmann  Equation 

The  Boltzmann  equation  of  a  single  component  gas  is  given 
by: 

(2.1)  (■£-   +  i--^-)    f  =  J/   (f*f '"  ff»)B(0,U)de  d0  d^ 

where 

(2.2)  U  =  ||*  -  £|  . 

The  notation  is  that  usually  found   '    .   f(x,|,t)   is  the  mass 
density  in   (x,£,t)   space,   £  =  (i^i^A-*)      denotes  the  molecular 
velocity  and   £#   the  velocity  of  a  particle   after  a  collision.' 
A  star  subscript  is  used  instead  of   1  ,   in  order  to  avoid 
confusion  with   |,  .   The  function  B(6,U)   under  the  integral, 
is  representative  of  the  intermolecular  force,   with  9        and 
e   giving  the  orientation  of  a  collision. 

We  will  consider  a  gas  which  is  slightly  removed  from  a 
state  of  equilibrium,  given  by 


-^2/2RT 
(27TRT  Y/d 


and  which  is  itself  a  solution  to  equation  (2.1).  p  and 
T  are,  respectively,  the  constant  density  and  temperature. 
A  small  departure  from  equilibrium  is  meant  in  the  sense  that: 

(2.4)    ff  =  f°(l  +  g) 
g2  «  1  . 


9.   H.  Grad,  "Principles  of  the  Kinetic  Theory  of  Gases," 
Handbuch  der  Physik,  vol.  12  (Springer,  1958) 
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For  a  problem  such  as  sound  propagation,  this  approxima- 
tion may  be  used  with  confidence.  Introducing  (2.4)  into  the 
Boltzmann  equation,  we  get  on  linearization 

(2.5)    (|t+  i-   !^)g  =  ^  /  f°(g'  +  g;  -  g  -  g*)B(e,u)de  de  ai. 

We  can  show  some  general  properties  of  the  integral  operator 
of  equation  (2.5)  by  introducing, 


(2.6)  j((J))  =   f  fuf"  [<\)}  B(e,u)de  de  d^  , 


with 


(2.7)        [(j)]  =  <i>;  +  r  -  <|>#  -  4>  . 

In  this  notation  it  can  be  shown 

(2.8)      /  if   J((|))d£  =  -^  /  f°f°[^][(t)]B  de  d0  df-*d£  . 

We  introduce  the  eigenf unctions,  § ,    of  the  linear  operator 
J,  by  writing 


o, 


(2.9)      J((J))  =  Af 


The  Maxwellian   f   is  extracted  from  the  eigenfunction  mostly 

7 
as  a  matter  of  convenience.   As  Wang  Chang  and  Uhlenbeck  have 

shown,  the  eigenvalues  are  real  and  non-positive.   For  from 

equations  (2.8)  and  (2.9)  we  have, 


10.   This  may  be  found  in  Ref.  1,  p.  67  and  in  Ref.  12,  p.  237< 
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,  r   f°f°\[fy]\2   B(0,U)d£  d9  d£*  d^ 

(2.10)    A=-iJ i- 5 

4     /  f°\$\2   d£ 


and  since  B(e,U)   is  positive   the  eigenvalue  A  is  non-positive, 
If  (J)  is  either   1,   4,   or   ^  ,   A  is  zero.   This  follows 
from  the  conservation  of  momentum,  energy  and  the  number  of 
particles,  in  a  two  "body  collision. 

Further  from  (2.8)  we  get  the  orthognality  of  eigen- 
functions  with  respect  to  the  weight  function  f   .   For  eigen- 
values  (j.,v   such  that   p,  ^  v  ,   we  have 

(2.11)  Cm-  -  v)  /  f°  $  f  di,  =   0 

where     and  if     are  the  corresponding  eigenf unctions.   If  we 
choose  to  represent  the  linear  operator   J ( (f) )   as 

(2.12)  J(<{))  =  /  KU,e*)<I>U*)de* 

we  immediately  see  from  (2.8)   that   K  is  a  symmetric  kernel. 
The  latter  property  is  usually  shown  by  construction  '  . 

Without  loss  of  generality,  we  may  take  the  eigenfunctions 
at   J  to  be  normalized  to  unity,  and  we  further  presume  that 
the  eigenfunctions,  which  we  denote  by   )(J).  I  ,   form  a 
complete  orthonormal  set.   We  expand  the  perturbed  distri- 


11.   Ref.  1,  p.  129 
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tuition  function  in   this    set, 


(2.13)  g  =  21     a     <|) 

n=o 


where  the  coefficients  are  given  by 


(2.14)    an  =  /  f°  g(()n  d| 


where   the  bar  refers  to  the  complex  conjugate.   As  we  will 
see  later  there  seem  to  be  several  "sensible"  ways  in  which 
to  order  the  eigenf unctions .   We  will  for  the  moment  think  of 
them  as  being  ordered  in  relation  to- the  ascending  order  of 
the  magnitudes  of  their  corresponding  eigenvalues.   If  the 
expansion  (2.13)  is  substituted  into  the  linearized  equation 
(2.5)  we  have 


(2.15)     (■*!+£•  -5^)  51   a  $     =  y      a   A  <l> 

v    -" '  K  at       ox7  * n  rn       *■ — p   n  n  rn 

n=o  n=5 


where  we  have  denoted  the  eigenvalues  by   A   .   One  should 
notice  that  the  summation  of  the  right  hand  side  starts  with 
n  =  5?  since,  as  already  mentioned  above,  there  are  five  zero 
eigenvalues.   By  making  use  of  the  orthogonality  property  we 
may  represent  equation  (2.15)  by 

(2.16)    -^  a.  +  b. .  *  3  a.  -  A. .  a.  =  0 

v     '  dt   1    ij   dx   j    ij   j 

where 
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(2.17)  b±J  =  b.±   =   J   f0^.  4  d^ 

and 

(2.18)  A.  .  =  (A.)5.  . 

The  parenthesis  In  the  last  relation  signifies  that  the  sum- 
mation convention,,  which  is  used  everywhere  else,  is  not  to  "be 

12 
used.   Since  the  matrix   b.  .   ,  as  Is  shown  by  (2.1?)  is 

hermitian,  the  system  of  partial  differential  equations  (2.16) 

is  of  the  symmetric  hyperbolic  type  discussed  in  the  litera- 
ls 
ture   „   The  system,  however,  has  an  infinite  number  of  variables 

and  equations.   In  later  sections  use  of  the  hyperbolicity 

will  be  made,  however  it  will  not  extensively  be  applied  in 

•this  study.   A  further  remark  in  this  connection  is  worth 

making.   The  "energy"  which  is  usually  defined  when  talking  of 

symmetric  hyperbolic  equations  is  the  sum, 


a. a. 

i  i 


In  this  case  however  this  turns  out  to  be  the  linearized 

14 
H-function 

We  now  leave  these  very  general  remarks  and  for  the  moment 


12.  b.  .  should  be  thought  of  as  a  three-component  vector,  each 

1J  component  of  which  has  the  two  scripts   ij. 

13.  K.  0.  Priedrichs,  Comm.  Pure  Appl.  Math.   J_  (1954);  H 
(1958). 

14.  This  has  already  been  pointed  out  by  K.  Grad,  ref.  3,  p. 
109. 
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■i  o  h 

at  least  take  B  =  B(0).   It  is  known  '    '   that  Maxwell  mole- 
cules (i.e.,  molecules  which  behave  with  an  intermolecular 
force  law  which  is  proportional  to  the  inverse  fifth  power  of 
the  distance  between  particles)  behave  according  to  such 
relation.   As  we  will  see,  this  assumption  greatly  simplifies 
the  analysis. 

For  Maxwell  molecules  then,  equation  (2.5)  becomes 

(2-19)     (Ir+ft-  |j)g-L(g) 

=  1  /  f°(g'+  g*  -  g  -  gjB(e)de  de  d^  . 

The  eigenfunctions  of  L(g)  are  the  complete  set  of  func- 
tions given  by  the  product  of  Sonine  polynomials  with  spherical 

harmonics.   Although  the  facility  offered  by  these  functions 

15 
was  noticed  for  instance  by  Burnett   ,  the  eigenfunction 

property  was  first  pointed  out  in  general  by  Wang  Chang  and 

Uhlenbeck  in  their  pioneering  report. 

Modulo  normalization,  the  eigenfunctions  are  represented 

by 


15.  D.  Burnett,  Proc .  Lond.  Math.  Soc .  40,  382  (1935) • 

16.  Wang  Chang  and  Uhlenbeck  show  this  in  Ref.  6.   It  was 
shown  independently,  although  somewhat  later,  by  L. 
Waldmann,  "Transporterscheinungen  von  mittleren  Druck", 
Handbuch  der  Physik,  vol.  12  (Springer,  1958)-   An 
elegant  and  simple  proof  of  this  property  is  given  by 
H.  Grad  in  Ref.  9. 
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where  the  molecular  velocity  is  normalized  with  respect  to  the 
equilibrium  thermal  speed, 

(2.21)     ?  =  ^/^/~RT^   . 

Referring  once  again  to  a  general  molecular  force,  Wang  Chang 
and  Uhlenbeck'  showed  that  the  eigenf unctions  of  equation  (2.6) 
always  have  the  angular  dependance  in  velocity  space  given  "by 
(2.20).   One  should  note  that  our  normalization  of  the  velocity 

o 

is  slightly  different  than  that  chosen  "by  Mott-Smith  ,  and 
Wang  Chang  and  Uhlenbeck'  (we  use  vART   and  they    J"2RT  ) . 
This  results  in  slight  variations  in  the  formulae  which  are 
given  on  the  next  few  pages.   The  Sonine  polynomials  can  be 
represented  by 


(2.22)     Sn(x)=-^i^  -^  /(l-s)m+neSX  \ 

ds11   - 


s=0 


For  a  discussion  of  a  Sonine  polynomials,  and  Legendre  polynomi- 
als  P*1  ,   one  is  referred  to  one  of  the  standard  works 
The  eigenfunctions  have  the  orthogonality  property  given  by 

^2  /s 

(2-25)    Y—  I   e_l  /2  Wr'i'm'^  =  5rr'&^'5mm.Nr^m 

(2tt)3/2 

where, 


17.  For  instance,  W.  Magnus  and  F.  Oberhettiger,  Formulas 
and  Theorems  for  the  Special  Functions  of  Mathematical 
Physics"  ,  (Chelsea,  194b). 
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£+1 

(o   ok)  N    =  2    P(r+i+3/2)  .  (J+\m\)l 

v  '  '  x  rim    — 

^77  '  r.'(2i+i)      (i-|m|  )J 


where  the  bar  again  refers  to  the  complex  conjugate.   Of  great 
aid  in  the  succeeding  work  is  the  following  recursion  relation 

(2'25)     (21+1)1-^  =  (i-m+l)[*rji+1,m-  4>r_i+1,m] 

+  2(i-»n)[(r+i+|)^i„lfm-  (r+l)(t)r+1>i_ljm]   . 

On  expanding  the  perturbed  distribution  function  in  terms 
of  the  eigenf  unctions ,    (j)  «   ,    we  get, 


(2.26)     g  =  21  a  .   * 


rim  yr£m 


with  the  relation  for  the  coefficients  given  by 


(2"  27)     ar£m-  pl~ -;  f°g  ^  d* 


Note  that  all  coefficients  are  immediately  dimensionless. 
Next  denoting  the  eigenvalues  of  L  by   ^p™      >      an^ 
substituting  the  expansion  into  the  linearized  equation  (2.6), 
we  get, 


-  17  - 


L(L*m      :  Kt>m       t> 


vim 


rim   Yrim 


(2.28) 


Kim    =   S 


rinr  rim 


rim   rinr  rim 


7 
Wang  Chang  and  Uhlenbeck  ,  have  shown  that   A  „   ,  regardless 

rim     ° 

of  the  molecular  model,  Is  Independent  of  the  index   m 
Therefore  each  eigenvalue  is  at  least (2i+l ) -fold  degenerate.   A 
further  degeneracy  for  Maxwell  molecules  is  pointed  out  in  the 
following  section.   The  latter  equation  has  been  made  dimension 
less  by  means  of  a  frequency   v  ,   in  the  following  way, 


f 


vt  =  t 


X1 


=  X 


(2.29) 


{        W, 


A 


rim 


=  A 


rim 


V_ 


We  will  leave  the  frequency   v   unspecified  for  the  present. 

Our  considerations  will  now  be  restricted  to  one-dimension- 
al unsteady  flows.   This  will  change  the  convective  terms  to 
£-,  ~^r   .   The  variable  x  now  represents  a  single  space  coordi- 
nate.  With  this  restriction  we  see  that  the  elements  of  the 
coefficient  matrix    b.  .  ,   defined  by  equation  (2.17)   is 
easily  evaluated  by  means  of  the  recursion  relation  (2.25). 
Examination  of  this  recursion  relation  shows  that  the 
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dependence  indicated  by  the  index  m  drops  identically  from 
equation  (2.28).   From  this  we  can  conclude  that  if  the 
initial  data  of  a  problem  is  homogeneous  with  respect  to  the 
index  m  ,   then  the  solution  will  not  develop  an  index  m 
in  time.   In  particular  if  the  initial  disturbance  is  isotropic 
in  the  plane  normal  to  the  x  direction,  we  can  set  m  =  0 
in  the  above.   This  latter  case  can  be  regarded  as  appropriate 
for  the  study  of  sound  propagation  and  will  be  the  case  which 
is  taken  up  here.   However  it  should  not  be  presumed  general 
enough  for  the  study  of  all  one-dimensional  problems.   For 
instance,  the  Rayleigh  problem  is  strictly  one  dimensional, 
but  nevertheless  requires  shear  moments  for  an  analysis. 
Henceforth  we  restrict  attention  to  the  m  =  0   case  and  denote 
the  eigenf unctions,  and  eigenvalues,  respectively,  by  (|)  ^ 

and   A  .  . 

Rather  than  use   (J)  »      we  will  with  Wang  Chang  and 
Uhlenbeck'  use  the  normalized  eigenf unctions, 

t  ^  ^ 

(2-30)    ifvl:   ■ 


2^1  (r+i+3/2) 
(tt1/2  ri(2i+l) 


18 
with  the  property 


18.    The  eigenfunctions  are  real  because  of  the  suppression 
of  the  index  m  . 
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(2. 3D  /   d|  co  yri*rli,    =    5rr,6^, 


with 


(2.^2)  03     =  -^77T 


(2TT) 


372 


The  recursion  relation  now  has  the  form, 


(2.33)  {2£+l)^,1ifr£  =   (£+1)     )frUl 


\ 


r£  +   3/2)2 


2i+l)(2i+3) 


•/ 


2r 


r-l,i+l 


(2i+3)(2i+l) 


+   iH 


ri-1 


W 


(r+£  +   1/2)2 
(2i+l)(2i-l) 


:- 


r+1    £-1 


2(r+l) 


(2£-l)(2i+l 


In   the   new  notation   the   perturbed   distribution   function   g   is 


(2.34)  g  =  ]E  *ri   *ri 


with 

(2.35)  b    »   =   /  a)  g  tf    .   d£ 


ri 


r^ 


And  the  linearized  equation  (2.15)  now  becomes, 


(2-56)    ^  +  h  ^r)  21  *ri*ri 


bri  Ari^rJ 
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Several  examples  of  the  eigenfunctions  and  their 
associated  coefficients  are: 


(2.37) 


r 


if  _  =  1 


oo 


*01   =    *1 


< 


tf 


10 


> 


i 


02 


4   (i-^) 

2       v  J>J 


-G\e?  -  £ 


2 


'1 


3 


* 


11 


^ 


fV1 


3 


b01   =   Ul 


10 


sfP 


02 


11 


-      P 


11 


The  quantities  on  the  right  have  been  made  dimensionless 

with  respect  to  T  ,  JRT   ,   p   ,   and  are  respectively, 

o  ^   o     o 

the  perturbed  density,  velocity,  temperature,  stress  and 
heat  conduction.   In  the  present  notation  the  "gas  law" 
identity  is 

(2.38)    p  -  p  +  T  . 

Here   p  is  of  course  the  dimensionless  pressure,  and  it  will 
be  eliminated  in  what  follows  by  means  of  relation  (2.38). 

We  now  wish  to  generate  differential  equations  for  the 
determination  of  the  coefficients,   b  »  .   First  substituting 
the  recursion  relation  (2-33)  into  equation  (2.36)  we  get, 


21  - 


(2.39)  0 


=  l^^ItT   -   \*Kt 


+ 


ri 


(l+i)  <  i 


v,H+l 


2l+T7(2l+J7 


-  if. 


2v 


r-l,i+l    /(2i+3) (2£+l 


+ 


,r+l+    1/2  )2~^ 
V^-l    /(2i+l)(2i-i; 


2(r+l)       ' 
Vr+i,.£-l    V (2i-l)(2i+l) 


And  finally  operating  on  the  above  with 


(2.^0)     /  di   co  i 


U.V 


we   get   the    infinite    system  of  differential   equations   governing 
the   coefficients,    b 


(2.41)  0    =    (^rr   -    A'   ,)b    „   +   | 

v  '  Mt  '  |_iv '    u.v        d 


:■: 


,,[i+v+    1/2)2    ' 
bH,v-l  yt2v-l)(2vll) 


-    b 


2\i+2 
\l+1,v-1   )/  (2v-l)  (2v+l 


HwlLipS 


-   b 


2p. 
p.-l,V+l   V  (2v+3)  (2v+l 


As   will    be   pointed   out    in   the   next    section,    the    smallest 
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non-zero  eigenvalue,  in  magnitude  (all  non-zero  eigenvalues 
are  <  0,  see  (2.10))  is   A-,-,   and  we  normalize  with  respect 
to  it,  i.e.,  we  take 


(2.42) 


11 


in  the  normalization  given  in  (2.29).   Using  the  relations 
given  in  (2.38)  and  the  fact  that  the  first  three  eigenvalues 
vanish,  the  first  three  equations  of  (2.4l)  are 


<2-^)   ^rrP+sn-«  =  ° 


2.W)   ^f  +  ^r  P1X  +  £r  (P  +  *) 


=  0 


to     h^\  ^  ,  2     S   q      2   S   n 


0 


These  will  be  recognised  as  the  linearized  and  dimension- 
less  conservation  of  mass,  momentum,  and  energy  equations. 

As  is  well-known  this  system  is  not  determinate,  it 
contains  more  unknowns  than  equations.   Which  is  in  fact  the 
case  for  any  finite  collection  of  the  equations  of  (2.4l). 
If  the  system  given  in  (2.43-2.45)  is  terminated  by 
excluding  the  quantities  which  do  not  have  a  time  derive  we 
have,  on  returning  to  the   "b"  notation, 


(2.46) 


1 

0 

0 

0 

1 

0 

0 

0 

1 

&♦ 


0         1 

0 

1         0 

>2/3 

0    -^2/3 

0 

—               — 

b 
00 

b01 

bio 

l_         L 

0 
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These  are  just  the  inviscid  equations  of  gas  dynamics,  but 
arranged  such  that  the  coefficient  matrix  of  the  space 
derivative  is  symmetric.   This  as  was  shown  earlier  is  a 
general  property  of  any  truncated  system  derived  from  the 
full  set  of  (2.41). 

It  is  clear  that  the  ordering  of  the  coefficients   b  „ 
is  arbitrary.   A  complete  solution  of  (2. 41)  will  naturally 
be  independent  of  the  ordering,  but  it  goes  without  saying 
that  a  complete  solution  is  not  feasible.   On  the  other  hand 
an  approximate  solution,  for  instance,  one  gotten  by  the 
truncation  (2.4l)  will  often  be  strongly  dependent  on  the 
ordering  of  the   b  «    .      In  fact  various  orderings  are  appropri- 
ate for  different  phenomena.   This  will  become  more  evident 
as  we  proceed.   We  will,  however,  now  indicate  two  orderings 
which  seem  natural  and  as  we  will  see,  are  appropriate  to 
different  situations. 

The  first  ordering  is  that  which  is  almost  universally 
chosen,  according  to  the  degree  of  the  polynomial  f    »    .      We 
see  from  the  expression  given  in  (2.20)  that  f    „      is  a  poly- 
nomial of  degree   2r+i.   The  ordering  of   b  .   can  thus  be 
chosen  to  correspond  to  the  degree  of  f    *      and  subordered 
in  ascending  values  of  the  index  SL    .      For  example 

(2.47)    bQ0,  bQ1,  b1Q,  bQ2,  blx,  b^,  b2Q,  b12,  bQV  ... 

The  other  ordering  which  seems  natural,  is  in  correspondence 

-  24  - 


with  the  ascending  magnitudes  of  the  eigenvalues,   I  A  J, 
as  was  introduced  at  the  outset  of  this  section. 

Before  leaving  this  section  we  one  again  take  up  the 
case  of  non-Maxwell  molecules.   Aside  from  Maxwell  molecules, 
no  molecular  force  law  permits  the  explicit  representation  of 
its  corresponding  eigenf unctions.   This  is  really  less  of  a 
liability  than  it  seems.   In  the  next  section  we  consider  the 
spectrum  of  the  entire  linear  Boltzmann  operator.   As  will  be 
seen  we  shall  be  succesful  in  obtaining  certain  power  series 
representations  of  the  spectrum.   This  power  series  develop- 
ment is  crucially  dependent  upon  the  Maxwell  molecule  eigen- 
functions  exhibited.   More  specifically  It  depends  on  the 
recursion  relation  given  by  (2.25).   This  fortifies  the 
natural  tendency  to  expand  the  distribution  function  in  terms 
of  Sonine  polynomials  in  spite  of  the  fact  that  they  are  not 
in  general  eigenf unctions  of  the  integral  operator.   In 
regards  to  the  differential  operator  of  the  Boltzmann  e- 
quation  they  are  so  to  speak  the  natural  functions  to  expand 
in. 

The  one  dimensional  form  of  equation  (2.5)   will  now 
be  considered.   The  previous  argument  given  to  disregard 
the  shear  moments  still  aplies.   We  therefore  use  the 
expansion  given  by  (2.34).   This  gives, 
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<2-48)       <^+*i'ST>2l*ri*rJ 


y  t>ri  /c%l> i]B(0,U)de  d0  d?* 


where  B(6,U)   has  been  normalized  accordingly.   On  operating 
on  (2.48)  with  (2.40)  in  order  to  generate  moment  equations 
we  obtain  the  same  differential  operator  as  is  found  in 
equation  (2.4l).   Instead  of  a  single  term  appearing  from 
the  collision  operator  we  now  get, 

(2.49)  5~A        b        =  Y"  b         fco„ao[TA      ]f     B(0,U)de    dd    d?#d£    . 

v  '  £—  ru.v   rv        ^—     rv   J    *        rv      \iv 

The  contraction  of  the  second  subscript  follows  from  the 
remark  following  equation  (2.21).   For  Maxwell  molecules  the 

o 

A     vanishes  for  r  /  n  .   Mott-Smith  has  considered 
A     for  a  wide  range  of  molecular  force  laws,  and  has 
shown  that   A     is  "almost  independent  of  the  force  law, 
at  least  for  small  r,|_i,v".  It  is  in  fact  the  case  that 

o 

several  of  these  elements  are  independent  of  the  force  . 
This  being  the  case  we  can,  for  non-Maxwell  molecules  take 
with  some  confidence  only  the  diagonal  part  of   A  _,  . 
One  may  easily  show  that  this  furnishes  what  Chapman  and 
Cowling  call  the  first  approximation  in  the  determination 
of  the  transport  coefficients. 
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3-   Eigentheory 

This  section  will  be  devoted  to  the  computation  of  the 
eigenvalues  of  the  linearized  Boltzmann  operator;  and  the 
subsequent  calculation  of  the  eigenvalues  of  the  differential - 
integro  operator  (2.5).   The  latter  operation  is  equivalent 
to  the  determination  of  the  spectrum  as  given  by  the  dispersion 
relation. 

Prom  the  previous  section  we  may  represent  the  eigenvalues 
of  the  integral  operator  by 

(3-D  \£    -I   ^    co^L(^)  . 

All  but  the  0-integration  in  (3-1)  can  be  carried  out''   to  get 

(3-2)      7^  =  2tt  Ae  B(0)[cos2r+i  §  P^cos  §■) 

+  sin21^  I  P.  (sin  |)  -  (1  +  5   6  .)]    . 

2  £ v     2  ro  rl'- 

-3/2  ! >7 
For  small  values  of  9,   B(S)   behaves  as  0  '  ,     so  we  see 

that  individual  terms  of  the  integrand  lead  to  divergent  integrals 

However,  the  combination  of  terms  erases  this  divergence,  as  can 

be  seen  by  inspection.   We  see  from  (3-2)  that 

(3-3)  A   =  \     _   _  . 

v    '  ro    r-1,  1 

This  degeneracy  will  play  an  important  roll  in  the  succeeding 

7 
work.   Wang  Chang  and  Uhlenbeck  have  evaluated  some  of  the  eigen- 

o 

values  and  later  Mott-Smith   evaluated  many  more.   A  compilation 
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of  this  may  be  found  in  reference  17-   Table  I  contains  the 

collection  of  values  of  A  '   =  A  ./Xn ,  ,   taken  from  reference  8. 

vl  rf  11* 

One  notices,  as  Mott -Smith  points  out,  that  in  addition  to  the 
above  degeneracy,  we  have  "h-,„   =   Ap   =  A   .   This  triple  degeneracy 
is,  of  course,  in  addition  to  the  triple  degeneracy  A   =  A   = 
A   =  0.   The  latter,  as  pointed  out  already,  is  a  consequence  of 
the  conservation  of  mass,  momentum,  and  energy  in  a  two  particle 
collision.   We  will  in  the  future  refer  to  these,  and  the  branches 
which  they  define,  as  the  hydrodynamic  roots.   The  reasons  for 
this  will  become  clear. 

In  order  to  study  the  character  of  plane  wave  solutions  to 
(2.4l),  we  assume  a  dependence 


(3.4) 


bri  =  hr£ 


at ' -ikx 


for  all  dependent  variables  where  the   b  ,   on  the  right  is  a 
constant.   This  corresponds  to  a  Laplace  transform  in  time  and 
Fourier  transform  in  space.   However,  since  we  will  not  for  the 
moment  be  interested  in  initial  or  boundary  value  problems,  we 
do  not  formally  perform  the  transforms.   Making  the  substitution 
of  (3-4)  in  (2.41) ,  we  get 


(3-5)  (a-*      )b      -   ik 


(U+v+l/2)2  ~* 
Vv-1    V(2v-l)(vu+D 


"V+l,v-l  \(2v-l 


M.+2 


)(2v+l) 


hv+i^HtH? 


2m. 


"V-l,v+l    V(2v+3)(2v+l) 


0 
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Using  the  ordering  by  eigenvalues  we  can  write  the  system  as 


(3.6) 


[al  -  ikB  -  A]V  =  0 


where   I   is  the  identity  matrix,  and 


(3-7) 


A  = 


00 


01 


A 


10 


02 


B   is  the  coefficient  matrix  of   ik  and  V  is  a  vector  whose 

components  are  given  by  (2.47)   The  remainder  of  this  section 

/  \  19 
will  be  spent  in  locating  the  roots   a(k)     of  the  determinant 

of  equation  (5.6),  or  equivalently  finding  the  eigenvalues  of 

the  matrix 


(3-8) 


C  =  +ikB  +  A 


The  determinant  of  the  linear  homogeneous  system  (3-6) 
will  be  referred  to  as  D(a,k).   This,  of  course,  defines  the 
dispersion  law.   To  avoid  questions  of  rigor  we  will  consider 
the  system  (3-6)  as  arbitrarily  large  but  finite.   Bu  this  dodge 
we  may  use  the  results  of  the  function  theory.   A  truncated 
matrix  of  the  system  (3-6)  is  given  for  the  eigenvalue  ordering 
in  Table  II  and  for  the  polynomial  ordering  in  Table  III. 


19-   For  a  boundary  value  problem,  one  is  essentially  interested 
in  k  =  k(a).   This  case  is  considered  in  the  reference  of 
footnote  7. 
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On  setting   k  =  0   in 

(3-9)   ■   D(k,a)  =  0 

we  see  immediately  that 
(3-10)     a   =   X^ 

are  the  roots.   As  already  noted  a  =  0   and   a  =  >>„p  (=  A   =  ^01 ) 

are  triple  roots  and   a  =  A  ..(  =  A     )   are  double  roots.   We 
^  rOv    r-1,1' 

will  now  develop  the  roots  in  a  power  series  in  k   in  the  neigh- 
borhood of  each  of  these   k  =  0  points.   Each  k  =  0  point 
defines  the  origin  of  one  or  more  of  the  branches  of  (3-9) •   That 

such  a  power  series  exists  follows  from  function  theory.   As  is 

2 
shown  in  Appendix  1,  the  dispersion  relation  is  a  function  of  k  . 

In  the  neighborhood  of  a  simple  zero  the  root   a,   may  ve  developed 

2 
in  a  series  of  ascending  powers  of   k  .   However,  in  the  neighbor- 
hood of  a  double  or  triple  point  a  Puiseux  series  must  be  used. 
In  our  case  this  amounts  to  a  power  series  expansion  in  k.   For 
a  discussion  of  these  remarks  one  is  referred  to  a  book  on  function 

4-V,  20 

theory. 

We  write  for  the  roots   a, 

(3-H)     o   =   o„   +   a.,  ik  +  c„k  +  ... 


and  for  D(k,a)   we  make  the  expansion, 


20.   E.  Goursat ,  A  Course  in  Mathematical  Analysis,  vol.  II,  pt .  1, 
(Dover,  1959~T 
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(3.12,    »-».,♦ -£♦££♦  ££♦...  ■ 

The  primes  denote  differentiation  and  the  zero  subscript  the 
evaluation  at  k  =  0.   Relation  (3-ll)  is  placed  in  (3-12)  and 
this  results  in  the  evaluation  of  the   a   in  (3«ll)>  when 
(3.12)  is  placed  equal  to  zero.   The  latter  states  that 


(3.15)    o  -  »„  -  l£>-  1^-  ...   -i£>-  . 


It  is  shown  in  Appendix  1  that  for  our  particular  determinant 


(3.1>) 


dnD 


dk 


n 


D(a! 


hh 


1  j 

n^n 


1^! 


a.  ,  , 
±2j2 


k=0 


>       &\         t        ) 


1  J 
ndn 


k=0 


1 

As  explained  in  the  appendix,  the  symbol   D(a.  .)   represents 
the  determinant  with  its  (ij)  element  replaced  by  its  derivative 
In  general,  (3-1^)  represents  the  determinant  with  n   of  its 
elements  replaced  by  their  respective  derivatives.   And  when 
an  element  appears   m  times  it  represents  the  m   derivative  . 
of  that  element.   The  relation  given  by  (3.14)  will  be  a  great 
aid  in  the  study  of  the  dispersion  relation. 

In  the  succeeding  section  we  sill  treat  in  general  all 
branches  which  occur  in  the  expansion  of  D(c,k).   Worthy  of 
special  attantion  are  the  three  branches  which  occur  in  connec- 
tion with  the  o=0     triple  point,  the  hydrodynamical  case.   We 


31 


treat  this  next  and  in  greater  detail  than  any  other  branches 


The  Hydrodynamical  Branches 

In  this  case  it  is  more  convenient  to  examine  Table  III, 
the  polynomial  ordering,  rather  than  the  eigenvalue  ordering 
of  Table  II  which  will  be  used  in  all  other  cases.   (Of 
course,  the  results  are  indifferent  to  whatever  ordering  we 
use.   We  use  a  particular  ordering  only  for  the  convenience 
it  affords  us).   To  facilitate  remarking  on  Table  III  we  have 
placed  "abcissa"  and  "ordinate"  lines  marked  1,  2,  and  3  in 
the  table.   Further,  we  shall  refer  to  abcissa  as   i   and 
ordinate  as   j.   For  instance   1  <  i,j  _  2   refers  in  the 
table  to  the  square 


(3-15) 


a    -   AQ2 

/&"< 

/5- 

The  hydrodynamic  branches  are  found  by  taking 


(3.16)  o   =  ictk+  pk2  +   iyk-5  +    5k     +   iek5  + 


in   D(k,c),  i.e.,  in  Table  III.   With  the  substitution  of  (3-l6) 

d      d 
it  is  clear  that   D,  -57-  D,   — ~  D  vanish  identically  when  k 


■  1; 


dk< 


is  set  to  zero.   Hence  no  information  can  be  derived  until  we 
reach  the  third  derivative.   This  is  also  clear  from  the  fact 
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that  a  -  0  is  triple  root  at  k  =  0.   Our  formula  {~5.lk)    tells 
us  that  in  order  to  take  the  third  derivative  of  D  at   k  =  0, 
we  need  only  consider  determinants  in  which  three  derivatives 
occur  in  any  combination  (i.e.,  three  elements  differentiated 
once,  one  element  differentiated  thrice,  etc.).   Unless  the 
first  three  rows  (columns)  contain  a  differentiated  element 
the  row  (column)  vanishes,  hence  the  third  derivative  at 
k  -   0  at   D   is  equivalent  to  the  third  derivative  of  the 
3x3  determinant   i,j  <  1.   (in  this  case,  at   k  =  0   the 
remainder  of  the  determinant  is  composed  just  of  the  diagonal 
elements  which  are  eigenvalues   A.)   In  other  words,  for  the 
hydrodynamical  branch  the  dispersion  relation  is  given  by  the 
coefficient  of  k   in  the  expansion  of  the  first  3x3  deter- 
minant; we  denote  this  by 


(3-17)    coef  (ir5) 


o  -ik       0 

-ik     u     Y2jy  ik 
0   /2/J   ik     0 


0 


Computation  shows  that  this  leads  to 

(3.18)    a  -  +  fflS,    0  . 

One  will  observe  that  the  first  two  roots  are,  in  our  normaliza- 
tion, the  two  adiabatic  sound  speeds.   More  significantly,  one 
should  observe  that  the  determinant  given  in  (3-17)  is  precisely 
the  determinant  that  one  would  get,  in  our  normalization,  if 
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the  inviscid  gas  dynamic  equations  are  Fourier  and  Laplace 
transformed.   Hence  from  our  present  point  of  view,  that  of 
kinetic  theory,  the  inviscid  equations  of  fluid  dynamics  cor- 
rectly gives  the  dispersion  relation  to   0(k)   for  sound 
propagation  for  the  branches  which  we,  in  anticipation,  called 
hydrodynamical . 

The  above  paragraph  contains  the  reasoning,  in  detail, 
for  the  determination  of  a .   In  the  paragraph  below  we  will 
give  the  reasoning  in  somewhat  less  detail,  for  the  determina- 
tion of  p,  i.e.,  the  hydrodynamical  dispersion  correct  to 

4  2 

0(k  ),  and  the  roots  to   0(k  ).   After  this  just  the  results 

will  be  quoted  and  commented  on,  since  the  reasoning,  although 

not  complicated,  is  cumbersome  to  explain.   It  is  best  carried 

out  in  the  privacy  of  one's  thoughts  and  with  the  aid  of  (3-l4) 

is  quite  straightforward. 

4 
To  consider  the  determinant  to   0(k  )   it  is  convenient 

to  first  consider  the  possible  contribution  from  the  elements 

in  the  square   I ,  j  >  1 .   At  most  one  differentiated  element  can 

occur  in  this  square  (more  than  one  gives  zero  identically  for 

k  =  0) .   But  this  element  must  multiply  the  first  3x3 

determinant  differentiated  thrice,  which  vanishes  by  the 

results  of  the  calculation  to  0(k  ).   Therefore  the  calculation 

of  p   is  given  by 
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(3.19)    0  =  coef(k  ) 


a 

-ik. 

0 

0 

0 

-ik 

a 

Y2/y  ik 

-2//T 

ik          0 

0 

/2/3vik 

a 

0 

-/57f±k 

0 

-2/yy  ik 

0 

~A02 

0 

0 

0 

-Y5/-f  ik 

0 

-An 

For  a     we  must  stubsitute  each  of  the  three  possible  branches, 


(3.20) 


0      +  Pk^ 

l/5/3V  k  +  pk2 

-i/BT^  k  +  £k2 


each  of  which  in  principle  gives  a  different  determination  of 
p.   In  Appendix  1  it  is  shown  that   D(k,a)   has  only  real 
coefficients.   This  shows  that 


(3.21) 


D(k,u)  -  D(k,a*)  , 


i.e.,  if  a   is  a  root  then  so  also  is  its  conjugate  a*.   Hence 
ti  Y5/3     lead  to  the  same  p. 

It  is  important  to  note  that  the  determinant  in  (3-19) 
is  precisely  the  expression  one  would  get  on  transforming  the 
Navier-Stokes  equations  under  our  normalization.   The  last  row 
is  the  transform  of  the  relation  between  heat  conduction  and 
the  gradient  of  temperature;  and  the  fourth  row  the  relation 
between  stress  and  the  space  derivatives  of  velocity.   The 
'viscosity  and  heat  conductivity  coefficients  are  given  respec- 
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tively  by  the  reciprocals  of  A    and   A    (to  within  numerical 
factors),  in  our  normalization.   Again  from  the  point  of  view 
of  kinetic  theory,  the  Navier-Stokes  equations  give  the  hydro- 
dynamical  branches  of  the  sound  dispersion  relation  accurately 
to   0(k  ),   and  the  roots  to   0(k  ). 

The  calculation  of   y  is  furnished  by 


(3-22)      0  =   coef(k5) 


o 


•ik 


0 


0 


0 


•ik  a  /2/S  ik      -2//S  ik  0 

0  /27T  ik  a  -Yb/ 

0  -2//3X  ik  0 

0  0  -X575"  ik     Y17T5  ik  A1X 


c-A02  /8/15"  ik 


which  is  found  in  the  same  way  as  the  last  two  determinations 

The  results  of  the  determinations  indicated  in  (3-19) 
and  (3-22)  are, 


(3.23) 
(3.24) 


a  =  ^ —  +  o(k  ) 
11 


a  = 


+ 


i  VW5  [k  +  k3( 


+ -i+- 


5A02A11    lOA^    i5A^2 


)] 


+  k2(^+5f-  )  +  0(k4) 

;  11    ;  02 


The  branch  as  given  by  (3-23)  leads  to  a  purely  decaying  mode 
and  is  connected  to  heat  conduction.   Relation  (3.24)  furnishes 
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o 

the  two  propagating  modes.   With  the  addition  of  the  0(k  ) 
term  the  propagation  acquires  decay,  and  the   0(k  )  term  leads 
to  an  increase  in  phase  velocity. 

The  determinant  of  (3-22)  is  exactly  the  one  which  would 

2 
be  found  from  the   thirteen  moments"  of  H.  Grad  ,  under  our 

normalization.   (Of  course,  in  our  case  of  one  dimensional 

unsteady  flow  the  number  of  moments  are  reduced  to  five.) 

Actually  instead  of  a   appearing  in  the  last  two  diagonal 

elements  we  could  have  put  in  its  expansion;  this  was  done 

in  order  to  make  the  identification  with  the  moments  equations 

of  Grad.   As  before  we  say  that  the  "13  moments"  equations 

accurately  give  hydrodynamical  branches  of  the  dispersion 

relation  to   0(k  ).   One  should  note  that  the  Grad  equations 

actually  give  two  other  branches  of  the  dispersion  relation. 

This  will  be  dealt  with  in  Section  H- . 

It  is  furthermore  the  case  that  (3-22)  furnishes  the 

same  dispersion  relation  as  would  be  given  by  the  Burnett 

equations.   This  follows  a  procedure  equivalent  to  the  Chapman- 

Enskog  theory,  and  which  in  our  case  is  given  by  row  operations 

on  the  determinant.   Briefly  this  procedure  goes  in  the 

following  way.   We  multiply  the  fourth  row  of  (3-22)  by  a/^0p 

and  add  it  to  itself;  and  the  fifth  row  by  a/A    and  add  it 

to  itself.   The  last  two  rows  of  (3.22)  then  becomes 


(3.25; 


0   z£lk(i+  cr_)         0 

T5  \)2 


0 


0     -)/|ik(i  + 


-v 

yitik 

All' 

/gl* 

-\i 
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One  notices  the  absence  of  the  quadratic  terms  in  the  fourth 
and  fifth  columns.   These  do  not  appear  since,  as  pointed  out 
in  the  previous  paragraph,  only  terms  of  0(k)   enter  as  the 
coef  (k^)  in  the  expansion  (3.22).   The  remainder  of  the  reduc- 
tion is  carried  out  in  the  same  way.   The  mixed  terms   ka   are 
eliminated  by  multiplying  the  second  and  third  rows  of  (3-22) 
by  ik  and  suitable  constants  and  adding  to  the  fourth  and 
fifth  rows  as  given  by  (3-25).   Again  nothing  new  appears  in 
the  fourth  and  fifth  columns  by  the  previous  remark.   Lastly, 
the  (4,5)  and  (5,4)  elements  are  eliminated  by  the  suitable 
addition  of  row  5  to  row  4  and  vice  versa.   The  result  is, 


(3-26) 


k 


2 


0 


-2 


1/2 


02  ;   02         '     11 


k 


% —  +  4k — Jk     -  /5/3  ik 


JK 


11 


-y\ 


02 


02 


0 


0 


-A 


11 


which  are  the  transforms  of  the  Burnett  relations,  in  our 
normalization . 

Referring  to  Table  III,  and  the  determinant  given  by 

6 
i>j  <  3,   the  coefficient  of  k   leads  to   5,   and  the  coef- 

ficient  of  k   to   e   of  the  expansion  for  a      in  (3-l6).   The 

determinant  is  identical  to  the  one  which  would  be  gotten  from 

what  Grad  would  call  the  "26  moments"  equations.   In  this  way  we 


21.   C.  S.  Wang  Chang,  "On  the  Dispersion  of  Sound  in  Helium" 
APL-JHU,  CM-467,  UMH-3-F  (1948). 
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may  go  on  and  accurately  determine  the  dispersion  relation  of 
the  hydrodynamic  branches  to  any  desired  order.   This  property 
of  being  able  to  find  a      to  any  order  by  the  proper  truncation 
of  T)(o  ,k)      is  peculiar  to  our  determinant  and  is  traceable  to 

the  recursion  relation  (2.20). 

7 
Wang  Chang  and  Uhlenbeck   also  noted  this  property  of 

D(a,k)   and  carried  out  the  calculation  of  hydrodynamic 
branches  to  within   i,j  <  3   of  Table  III.   They  call   i,j  <  1 
"Euler",   i,j  <  2   "Navier-Stokes" ,   i,j  <  3   "Burnett",  and 
note  that  at  each  of  these  steps  the  dispersion  relation  is 
improved  by  two  orders.   They  also  note  that  their  dispersion 
relation  is  one  order  more  accurate  than  would  be  gotten  from 
the  equations  corresponding  to  the  appellation.   We  have  seen 
in  the  above  why  this  is  the  case.   A  more  appropriate  sequence 
of  appellations  would  be   "Euler",  "13  moments"  and  "26 
moments" . 

Before  going  on  to  the  remainder  of  the  dispersion  law  we 
pause  to  make  some  comments  on  the  hydrodynamical  branches. 
Although  we  found  the  transforms  of  the  equations  of  the 
Chapman-Enskog  procedure  in  the  above  analysis,  it  should  not 
in  any  way  be  construed  as  a  verification  of  the  Chapman-Enskog 
procedure  itself.   This  is  clear  from  the  fact  that  the  inver- 
sions of  the  relations  do  not,  without  proper  initial  conditions, 
give  Chapman-Enskog  relations  for  stress  and  heat  conduction. 
We  will  not  deal  with  this  in  this  report.   Next  we  note  that 
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the  expansion  of  a      is  in  inverse  powers  of  A,  i.e.,  with 
terms  of  the  type   (k/A   ).   As  we  take  larger  determinants 
we  gain  larger  A  „   and  we  can  presumably  describe  higher 
frequency  phenomena.   Unfortunately,  the  calculations  become 
more  tedious  with  larger  determinants.   It  should  be  noted 
that  within  the  framework  set  up  the  expansions  converge, 
however  slowly.   Again  for  the  reasons  given  above,  this  should 
in  no  way  be  construed  as  showing  convergence  of  the  Chapman- 
Enskog  procedure. 

We  now  turn  to  the  remaining  branches  of  the  spectrum  and 
determine  their  behavior  in  the  neighborhood  of  o  =  A  .  (k  =  0) . 
As  we  will  see,  this  will  not  necessarily  mean  that  the  wave 
number  be  small.   It  need  only  be  small  when  compared  with 
the  eigenvalue   A  .,      and  since  this  becomes  increasingly 
large  the  results  can  be  applied  to  larger  and  larger  wave 
numbers.   In  no  case  will  we  carry  the  analysis  as  far  as  in  the 
hydrodynamical  case  j  the  latter  was  given  a  more  lengthy  analysis 
only  because  of  its  special  significance.   We  will  separately, 
in  turn,  deal  with  the   cases   when  at  k  =  0,  o   is  non- 
degenerate,  doubly  degenerate,  and  triply  degenerate  (the  hydro- 
dynamical  case  is  an  example  of  this).   In  referring  to  D(k,a) 
we  shall  henceforth  mean  Table  II  only. 


ho  -Q'jjc:.crate  Bra ■  . e  s 

Since,  as  is  shown  in  Appendix  1,   D  =  D(a,k  ),  we  may 
express  the  branch  of  a  non-degenerate  root  as 
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(3-27)         a  =  A^,  +  Pk2  +  5k4  +  .  .  .  . 

With  the  help  of  the  relation  {~5.lK)    it  is  not  difficult  to 

see  that  the  determination  of  B   is  gotten  from  the  coefficient 

2 
of  k   in  the  expansion  of  the  determinant  in  which   k  =  0 

except  in  the  column  and  row  containing  \lv«   This  is  shown 

in  Table  IV.   If  this  is  carried  we  get 


(3-28)    B  =  - 


(V+1)  \(2v+l)(2v+3 


3T  *  -  A 

|iv    M.,v+1 


+ 


_2}± 1       1  +  v2  i  2,^+  X?l/2\  1  + 


(2V+3M2V+1)  A^-  Vl,v+lJ      }(2v+lM2v-l)  A^-  *— ^ 


2(M-+1) 


(2v-l)(2v+l)  A^-  \+1}V_± 


If  one  inserts  into  (3-24)  those  values  which  are  given  in 

Table  I,    it  is  found  that   8    is  negative.   It  is  plausible 

2 
to  presume  that  the  coefficient  of  k   in  the  non-degenerate 

branches  are  all  negative.   We  will  comment  on  the  signifi- 
cance of  this  in  Section  H- .   The  next  order  may  be  gotten  by 
filling  in  the  neglected  entries  of  Table  IV  and  expanding 
the  determinant. 


Doubly  Degenerate  Branches 

From  (3-3)  we  know  that  the  doubly  degenerate  eigenvalues 

are  given  by  A   =  A   _  _  .   Referring  to  Table  II  a  case  in 
&      J    ro    r-1,1  & 
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point  is   A   =  A   .   Using  the  same  line  of  reasoning  as 
used  in  the  preceeding  cases  we  see  that  a      of 

c   =   A         (=   A      )    +   ictk  +   (3k        is   given   by   the   coefficient    of 

2 
k   in  the  expansion  of  determinant   1  <  i,  k  <  2   in  Table 

II.   Carrying  out  this  calculation,  we  find 

(3.29)  a  =  +  2/yy  . 

To  compute   (3   which,  as  has  been  pointed  out,  is  the  same 
for  both  values,  we  need  only  consider  the  coefficient  of  k 
in  the  determinant  in  which  k  =  0,  except  in  the  two  columns 
and  rows  which  contain  A   ,  ^pr,-      0n  carrying  out  this  calcu- 
lation we  get 

<5-5°)     P  "  "  5A^  -  6(A20-A21)  -  30(A20-A12,  "  30(A20-A02)   ' 

The  general  expressions  for  a   and  P   may  be  derived  by 
returning  to  the  transformed  system  given  by  equation  (3-5)- 
If  the  equations  are  reordered  for  the  purposes  of  the  calcula- 
tion, we  may  write  the  matrix  in  the  neighborhood  of  A   (  =  A 
as  given  in  Table  V.    The  coefficient   a   is  immediately  seen 
to  be 

(3-31)         «  =  ±  K2r7T  • 

P   is  determined  by  taking  the  coefficient  of  k   in  the 
expression  of  the  determinant  of  Table  V.    This  gives 
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(3-32)     (3 


2r+3 


+ 


2r+l 


+ 


2(2r+3) 


v  rO   r,ly     ^  rO   r-1,0'      v  ro   r-l,2y 


+ 


2(2r-2) 
15(\.0-\-2,2J 


As  before,  the  next  term  in  the  expansion,  7,  is  furnished  by 
filling  in  the  remaining  entries  in  Table  V  and  expanding  the 
resulting  determinant.  This  is  not  done  since  the  result  does 
not  seem  worth  the  tedious  labor  involved.  The  (3  '  s  which  can 
be  computed  from  Table  I  are  negative.  It  is  again  plausible 
to  presume  this  negativity  for  all  p  given  by  (3-28). 


Triply  Degenerate  Eigenvalues   (^no  =   ^ 


A. 


02     21     30 ' 
We  shall  only  state  the  results  of  this  calculation  since 
it  offers  nothing  new  in  the  way  of  computation. 


0  =   So  "  k' 


(3-33) 


4 


3A 


+ 


8 


02 


W^-hl*      '  5l\32-W 


+  o(k4) 


•  =  A   ±  f?  ik  -  k' 


7 


O  -1  p 

15(A02-A12j  +  15(A02-A22j 


3 


+  ^N^W  +  2^02"A31J  j 


+  0(k3) 


Again  one  may  see, on  using  the  values  of  Table  I,  that  each  of 

the  modes  are  dissipative.   No  general  relation  is  given  in  this 

case  since  no  general  triple  degeneracy  is  to  be  found.   In  a 
sense  this  last  case  is  not  a  triple  degeneracy  at  all,  at 
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least  not  in  the  same  way  as  the  hydrodynamical  branches, 
is.   One  can  easily  see  that  the  branch  corresponding  to 
7\  may  be  -computed  as  if  it  were  a  non-degenerate  case. 

Further,  the  branches  corresponding  to  ~^-?n>    ^P1   may  be 
computed  as  if  they  were  a  double  degeneracy. 
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4.   Some  Remarks  on  Solutions  of  the  Boltzmann  Equations . 

Using  the  data  of  the  last  section  we  have  plotted  In 
Figure  1  the  roots  of  the  dispersion  relation  In  the  complex 
a-plane.   For  these  curves,   k   occurs  as  a  parameter  and 
the  dependence  of  a     on  k  may  be  found  by  returning  to 
the  appropriate  calculation  in  Section  J>.      Our  work  so  far 
has  been  aimed  at  solving  an  initial  value  problem.   The 
purpose  of  this  section  is  to  avoid  actually  solving  an 
initial  value  problem.   We  will  do  this  mostly  by  studying 
Figure  1,  and  making  some  plausible  assumptions.   In  this 
way  we  will  be  able  to  display  the  character  of  a  solution 
to  the  Boltzmann  equation  as  it  evolves  in  time.   This  will 
all  be  subject  to  a  somewhat  serious  restriction.   Our 
representation  of  the  dispersion  relation  was  in  a  power 
series.   Therefore,  if  only  a  few  terms  of  the  series  are 
taken,  the  representation  will  presumably  be  valid  for  only 

small  k.   Hence  we  will  be  restricting  ourselves  to  problems 

22 

with  sufficiently  smooth  initial  data   or  equivalently  that 

the  observer  point  be  far  removed  from  the  initial  disturbance, 

Looking  at  Figure  1  we  see  that  three  branches  emanate 
from  the  origin.   In  the  last  section  we  referred  to  these 
as  the  hydrodynamic  branches.   From  (j>.2J>)   we  have  for  the 
branch  on  the  real  axis, 


22.   In  our  normalization,  "smooth"  has  the  meaning  that 
the  ratio  of  the  mean  free  path  to  the  wavelength 
of  the  phenomena  be  small. 


-  ^ 


(4.1)      a        -  (310k2   ,    P1Q  >  0  . 

This  dispersion  would  also  be  obtained  by  considering  the 
equation 

ft-2'    slu  =  pioiJu 

which  we  identify  as  the  diffusion  equation.   Therefore  the 
mode  given  by  (4.1)  is  strictly  diffusing.   Prom  (3-24)  we 
have  for  the  upper  branch 


(4.3)    a  =JVT  ik  •  p_  k2 


oo  oo 


This  mode  therefore  gives  a  propagation  with  the  adiabatic 


speed    I 5/3 '  (in  our  normalization)  and  a  dispersitive  part 
by  transforming  a  differential  equation  of  the  form 


iven  by  p   k     The  expression  (4.3)  could  have  been  gotten 
oo 


(4.4)    |j  +  r^Tji  |H  =  p      _a£  u 

v    ;      dt    >J  ^7  v   dx    Koo  -n  2 


Prom  this  it  is  clear  that  one  gets  a  "wave"  propagation  with 
the  speed   /B/7^   around  which  there  is  a  diffusion  governed 
by  p>    .Of  course  the  conjugate  hydrodynamic  branch  gives 

OO 

the  same  phenomena  only  propagating  in  the  opposite  direction. 

Each  of  the  other  branches  can  be  considered  in  the 
same  manner.   The  chief  difference  between  the  other 
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branches  and  the  hydro-dynamical  branch  lies  in  the  fact  that 

j 

they  each  start  off  with  a  negative  real  part.   For  instance, 
consider  the  second  cluster  of  Figure  1.   From  the  last 
section  this  is  given  by 

(4.5)     a  =  A2Q  +  -£-  ik  -  P2Qk2  +  ... 

P20  >  0  ,    A2Q  <  0  . 


Such  a  dispersion  relation  could  also  be  gotten  by  trans- 
forming the  following  differential  equation, 

r  h    c  \      du    2   d       ,       Q     3 
(4-6)      St    +  7=  3x"  u  "  A20U  =  P20  T" 2  U  * 

This  is  of  the  same  type  as  (4.4).   The  undifferentiated 

A2Qt 
term  may  be  eliminated  by  taking  u  =  v  e      ,   which  decays 

in  time.   This  would  result  in  the  same  equation  as  (4.4). 

Therefore  the  second  branch  leads  to  a  mode  which  propagates 

with  a  speed   2/J~~jP  while  diffusing  according  to  (3  „  , 

and  decaying  with  a  folding  time  given  by   A~   .   We  see 

from  Figure  1,  and  the  results  of  Section  3,  that  the  branches 

lead  to  an  infinite  number  of  modes  which  propagate  and 

23 
diffuse,  and  an  infinite  number  which  jiast  diffuse.   All  but 

the  hydrodynamical  modes  decay  with  a  time  faster  than  the 

collision  time. 

As  we  have  mentioned  in  the  beginning  of  this  section, 


23^ Wang  Chang  and  Uhlenbeck  (ref.  7,  p,  32)  were  aware  of  the 
possibility  of  more  than  one  sound  speed;  they  specifi- 
cally point  out  a  mode  corresponding  to  equation  (4.6). 
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we  were  not  successful  in  calculating  the  eigenvalues,  a, 
for  all  k  but  only  for  small  values  of  k.   Hence  our 
dispersion  relation  is  only  capable  of  dealing  with  problems 
which  have  relatively  smooth  initial  data,  or  equivalently 
for  arbitrary  initial  data  the  observer  point  must  be 
sufficiently  far  from  the  location  of  the  disturbance.   On 
the  other  hand  the  dispersion  law  is  capable  of  dealing 
with  problems  for  all  time.   If  one  wishes  to  consider 
phenomena  for  times  larger  than  the  collision  time  the 
hydrodynamic  branches  need  only  be  considered,  since  all 
other  modes  are  asymptotically  small.  (Here,  as  in  previous 
discussion,  we  exclude  the  presence  of  boundaries).   And 
so  for  sufficiently  large  times  and  smooth  initial  data 
(or  the  equivalent  given  above)  we  get  the  same  modes  as 
given  by  Navier-Stokes ,  Burnett,  etc.  (see  identifications 
made  in  Section  3)-   This  again  should  not  be  construed  as 
a  derivation  of  these  systems  of  equations,  since  we  have 
not  discussed  the  way  in  which  initial  data  enters. 

One  needs  not  only  consider  the  hydrodynamic  branches. 
The  modes  of  all  the  branches  to  the  right  of  some  arbitrary 
point  on  the  real  of  Figure  1  may  be  considered,  and  all 
the  modes  to  the  left  would  be   considered    small  when 
compared  to  the  modes  to  the  right.   Indeed  the  entire 
solution  including  all  the  modes,  may  be  thought  of  as  a  sue 
cession  of  epochs,  or  temporal  boundary  layers.   Each  epoch 
is  determined  by  the  "folding  times"   A",  .   After  each  of 
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these  times  a  cluster  of  modes  becomes,  roughly  speaking,  e~ 
of  its  initial  value.   This  does  not  at  all  mean  that  it  is 
negligible  when  compared  to  the  modes  immediately  to  the 
right  of  it  in  Figure  1.   Inspection  of  Table  I  shows  that 
ratios  of  eigenvalues  need  not  be  that  large  to  warrant 
such  a  statement.   However  a  particular  mode  can  be  negli- 
gible compared  to  modes  sufficiently  far  enough  to  the 
right  of  it  in  Figure  1  .   A  cluster  of  modes  is  washed  out 
of  the  picture  more  quickly  than  its  neighboring  cluster 
to  the  right  and  less  quickly  than  the  one  to  the  left.   The 
hydrodynamic  mode  is  the  last  of  these  epochs,  and  all  other 
modes  are  asymptotically  small  compared  to  it. 

Lastly  we  would  like  to  say  a  few  words  about  the 
infinitude  of  sound  propagations.   It  is  indeed  the  case 
that  the  infinity  of  sound  propagations  as  they  appear  in 
this  report  are  specifically  a  product  of  having  chosen 
Maxwell  molecules.   In  general  we  have  no  a  priori  reason  to 
presume  that  the  eigenvalue  degeneracy  of  Maxwell  molecules 
is  likely  to  occur  for  other  molecular  models.   However, 
there  are  several  arguments  which  can  be  offered  in  support 
of  the  extraordinary  sound  speeds  in  general.   We  recall  from 
Section  2  that  from  the  isotropy  of  the  Boltzmann  integral 
operator  Wang  Chang  and  Uhlenbeck'  were  able  to  show  that  a 
(2^+1 ) -fold  degeneracy  of  the  eigenvalues  occurs  for 
arbitrary  molecular  force  laws.   This  degeneracy  was  wiped 
out  in  our  problem  by  the  elimination  of  shearing  moments. 
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As  pointed  out  in  Section  2,  however,  a  general  attack  on  one 
dimensional  unsteady  problems  must  include  such  shearing 
moments.   Therefore  a  degeneracy  would  occur  and  we  could 
expect  propagation  arising  from  it. 

Still  another  argument  in  favor  of  the  infinitude  of 
speeds  is  furnished  by  the  remarks  made  at  the  close  of 
Section  2.   For  non-Maxwell  molecules  we  no  longer  have 
the  degeneracies  which  were  responsible  for  the  non-hydro- 
dynamic  propagations  found.   However  from  the  remarks  at 
the  end  of  Section  2,  we  know  that  for  other  molecular 
models  the  trace  elements  of  the  collision  matrix  operator 
are  only  slightly  removed  from  the  Maxwell  values.   The 
double  and  triple  points  are  replaced  by  adjacent  sets  of 
points.   In  the  next  section  we  will  show  how  these 
neighboring  points  can  lead  to  a  degeneracy  and  give  wave 
propagation.   The  argument  will  rest  on  the  fact  that  each 
truncated  system  of  moments  equations  is  hyperbolic.   This 
has  the  effect  of  producing  wave  propagation  for  sufficient 
ly  high  frequency  phenomena.   Analytically  what  can  happen 
is  that  the  adjacent  points  lead  to  diffusion  for  small  wave 
numbers  but  at  some  threshold  value  they  coalesce   and 
produce  propagation. 
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5.   Kinetic  and  Macroscopic  Models 

Throughout  this  report  we  have  concentrated  our  efforts 
in  obtaining  the  dispersion  relation  of  the  exact  equation, 
i.e.,  for  the  linearized  Boltzmann  equation.   These  lead  to 
series  representations  of  the  branches  of  the  dispersion  re- 
lation, as  given  by  the  expressions  found  in  Section  2. 
Such  expressions,  of  course,  limit  us  to  the  study  of 
relatively  small  wavelength  phenomena.   However,  in  principle, 
the  series  representations  may  be  carried  to  any  desired 
order,  but  unaided  by  machines,  such  calculations  soon  become 
discouraging.   In  the  hope  of  gaining  more  control  on  the 
entire  range  of  wavelengths,  we  now  give  up  the  notion  of 
finding  the  dispersion  relation  for  the  exact  equation.   We 
shall  rather  explore  the  possibility  of  finding  approximate 
formulations  which  lead  in  some  sense  to  approximate  dis- 
persion relations. 

Certainly  the  simplest  approximation  scheme  would  lie 
in  using  a  truncated  system  without  restriction.   For 
instance  referring  to  Table  III,  if  we  truncate  the  system 
i,j  <  1   off  from  the  infinite  system  we  have  the  following. 
There  are  three  roots  of  the  dispersion  relation.   Two 
correspond  to  propagation  in  opposite  directions  with  the 
adiabatic  speed,   1375  >      and  the  third  is  a  zero  speed 
propagation.   Equivalently  we  could  have  gotten  the  same 
results  by  truncating  the  system  given  by  the  differential 
equations  (2.4l).   (These  can  also  be  gotten  by  direct 
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inversion  of  the  three  by  three  matrix,   l,j  <  l).   The 
propagations  are  of  course  the  characteristic  speeds  of 
the  Euler  equations,  with  the  zero  speed  corresponding 
to  the  inviscid  limit  of  heat  conduction.   This  system 
is  free  of  dissipation  and  gives  a  propagation  speed 
independent  of  the  wave  number.   The  dispersion  relation 
is  given  parametrically  in  k  ,   in  the   a-plane  in  Figure 
2. 

The  next  truncation  which  is  implied  would  be  given 
by  i,j  <  2   in  Table  III.   More  directly  by  the  corresponding 
truncation  of  the  system  of  differential  equation  (2.4),  we 
get 

— 

3d    3 


ft  +  H=° 


a 


+  t|(p  +t>  +  Jpn  =  ° 


(5.1)^ 


2  oil    2  C  Sl    n 


St    +  3^  +   5~"3 


x 


apll   .         4  Su    8   aSl   n 
St~  +   A02pll  +"ox+T5+_5x""° 


^Sl    .   -   .  Spll  .  3    ST  _  0 

sr  -  Ansi  +  ~^r  +  2  ^  -  ° 


These  will  be  recognised  as  the  thirteen  moments  equations 

2 
of  H.  Grad  .   This  was  also  pointed  out  in  section  3,  when 
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considering  the  hydrodynamic  branches  of  the  spectrum.   In 
fact  in  the  method  of  truncation  we  have  been  tacitly  ass- 
uming that  we  are  following  the  truncations  demanded  by 
the  computation  of  the  hydrodynamic  branches  of  Section  5- 
A  shortcoming  to  this  method,  as  well  as  an  alternative 
plan  will  be  presented  shortly.   As  has  already  been  pointed 
out  the   "13-moments"   dispersion  relation  contains  the 
dispersion  relations  of  Navier-Stokes  and  Burnett. 

In  Section  2  it  was  shown  that  any  truncation  no 
matter  how  large,  gave  a  symmetric  hyperbolic  system  of 
equations.   This  has  several  consequences  which  are  of 
importance  to  us.   Rather  than  deal  with  the  "13-moments" 
equations  specifically,  we  treat  any  truncated  system  no 
matter  how  large.   First  we  mention  that  any  such  system 
has  real  finite  characteristic  speeds.   In  our  case  this 

means  that  as   k  -)  00  we  are  lead  to  finite  speeds  of 

propagation  (some  of  which  may  vanish).   Next  we  may 
estimate  the  dissipation  of  such  a  system.   We  represent 
a  truncated  system  by 


(5-2)      (I^  +  A^-A)b  = 


0 


where   I   is  the  identity  matrix,   A   the  coefficient 
matrix,  as  gotten  from  (2.4l),  and  /\     the  diagonal  matrix 
made  up  of  the  eigenvalues  of  the  Boltzmann  integral 
operator.   The  dispersion  relation  of  (5.2)  is  given  by, 
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(5-3)     det  [la  -  ikA  -  A]  =  0   , 

An  equivalent  problem  Is  to  find  the  eigenvalues   a  ,  such 
that 

(5-4)      (ikA  +  A)ha  =  aba 


where  b   is  the  corresponding  eigenvector.   Multiplying 

...  * 

(5.4)  by  the  conjugate   b   ,  b   we  get 


(5.5)     ik(A  ba  ,  b*)  +  (A  ba  ,  b*)  =  a(ba  ,  b*) 


Since   A   is  real  and  symmetric   (Ab,b  )   is  real.   Therefore 
if  we  denote  the  real  part  of   a  by   a   we  have 

*, 


(5.6)    Of  =  - 


<V  V 


which  is  a  very  revealing  expression.   Since     has  been 
shown  (see  Section  3)  to  be  non-positive,  we  have, 


(5-7)     ov   <  0  . 

* 

Further  if  we  denote  by   A   the  eigenvalue  of  largest 


magnitude,  we  have  further 


(5-8)     A*  <  ov   <  0 
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Hence  any  truncated  system  has  limited  decay,  and  as  the 
size  of  the  truncated  system  increases   A   becomes  more 
negative.   It  should  he  born  in  mind  that   a   is  a  function 
of  k,   as  is  the  eigenvector  b   .   The  work  of  the  previous 
sections  indicates  that   a   becomes  more  negative   with 
increasing  k.   A  statement  which  is  certainly  physically 
reasonable.   Further,  by  increasing  the  size  of  the  system 
it  is  plausible  to  assume  that   a  (k)   will  become  more 
negative.   This  is  clearly  implied  from  equation  (5.6).   We 
conclude  from  this  line  of  reasoning  that  a  truncated  system 
in  limiting  the  decay  of  high  wave  number  phenomena,  inade- 
quately describes  this  part  of  spectrum.   Also  in  support 
of  this  statement  are  the  calculations  of  section  3-   We 
showed  there  that  any  coefficient  in  the  series  representation 
of   a  is  not  correctly  given  until  a  large  enough  matrix 
is  taken. 

We  return  to  the  thirteen  moments  system  as  given  by 
(5.1).   The  dispersion  relation  is  given  by   i,j  <  2   of 
Table  III.   For  k  =  0  ,   the  roots  are   a  -  0,0,0,  ^11^02    • 
The  triple  root  at   (c,k)  =  (0,0)  has  already  been  analyzed 
for  small   k  ,   in  Section  3  and  leads  to 


(5.9) 


-  /57T  ik  +  k2(^-  +  ^-)   +  0(k5) 
;  11    -^  02 


c  =  ^  +  0(k4) 


11 
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Using  the  methods  of  Section  3,  we  immediately  get  for  the 
other  two  roots, 

,2  8A  k 

(5-iO)         c  -  A       -  £—  [5/3  +  1C-/A     ,A  -r  ]  +  0(k4) 

11         a1x  l3^AH      02 ; 


k         [il/3  02 , 

15 (*-,-, -'W 


(5-11)  V,9    -  f-  [4/3  +   ^/vlj    ^   -   °<k 

.1       02 


On   substituting  the   values   of   Table    I,    we   find   that   both 
branches   given  by    (5-10)    and    (5-H)    move   to  the   right.      This 
is   contrary  to  physical   experience    since   it   predicts   a 
decrease   in  decay  with  wave   number.      When 

(      k  ~  +    •    35 
(5-12)  < 

^     a  -   •    6AQ2 

these  two  branches  come  together  in  a  double  point  and 
leave  the  axis  in  a  way  shown  in  the  circles  region  of 
Figure  3.   Inspection  shows  the  disagreement  with  the  exact 
dispersion  relation  given  In  Figure  1.   This  is  directly  due 
to  the  polynomial  ordering  chosen.   We  shall  come  back  to 
this  point  after  sketching  In  the  remainder  of  the  "thirteen 
moments"  dispersion  relation.   In  order  to  examine  the  dis- 
persion relation  for  large   k,   we  extract   Ik   from  each 
row  of  the  determinant  given  by  (3-22J. 


-  56  - 


(5-13) 


X 

-1             0 

0 

0 

1 

x       )E73 

-2/^3 

0 

0 

/2?3          x 

0 

-/5/3 

0 

-2/yy    o 

x-e 

#1 

0 

0        -,[5X3 

s(8>15 

x-a 

0 


where 


(5-14) 


x  =    a/ik 

A02 

e   =      ik    ' 

ae   = 

ik 

One  then  considers   e   a  small  quantity  and  computes  the 
roots  in  this  limit.   For  Maxwell  molecules  using  the 
values  of  table  I,  the  case  which  is  sketched,  we  get 


(5-15)    o       +   5/9  AQ2  +  0(l/k^) 


(5-16) 


+ 


a   I  ik(^T5^)+  .33AQ2  +  0(l/k) 


+ 


-   ik(>/7olj")  +  -22A   +  0(l/k) 


The  coefficients  of   ik   in  (5. 15)  and  (5.16)  should  be 
recognized  as  the  "13  moments"  characteristic  speeds  . 
On  carrying  out  the  expansions  of  (5.36)  to  several  more 
terms  we  find  that  the  first  expression  is  to  be 
identified  with  the  hydrodynamic  branch  and  the  second 
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with  the  extraneous  sound.   This  completes  the  picture  given 
in  Figure  .3° 

We  pause  now  to  carry  on  our  argument  begun  in  Section 
h.      At  that  time  it  was  stated  that  although  a  system  may 
not  have  degenerate  eigenvalues  it  would  still  lead  to 
propagation.   This  being  because  the  truncated  systems  are 
always  hyperbolic.   The  mechanism  by  which  this  takes  place 
is  clear  from  Figure  3,   If  in  the  general  case  the  origin 
of  two  branches  are  close  together,  we  can  expect  them  to 
join  in  the  manner  shown  in  the  circle  of  Figure  3.   Instead 
of  moving  to  the  right  with  wave  number  they  will  probably 
move  to  the  left  however. 

We  recapitulate  the  remarks  made  in  connection  with 
the  thirteen  moments  equations.   For  smooth  phenomena, 
k  «  1  ,   they  accurately  describe  hydrodynamic  branch, 
but  not  so  for  the  branches  near   A-,-,  and   AQ?  .  Compari- 
son with  Figure  1  shows  that  the  dispersion  relation  is 
completely  distorted  in  that  region.   The  most  evident 
inconsistency  being  that  decay  decreases  with  wavenumber. 
For  large  wave  number  phenomena  the  inconsistency  is  more 
striking.   From  what  has  been  said,  and  also  some  remarks 
which  occur  later,  the  indication  is  that  all  branches  of 
the  dispersion  relation  move  toward  negative  infinity  with 
increasing  wave  number.   The  thirteen  moments  system,  or 
for  that  matter  any  truncated  moment  system,  has  a  finite 
absorption  cut-off.   Further  aggravating  this  is  the  fact 
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that  there  is  a  branch  for  which  the  absorption  decreases 
as  k  becomes  unbounded,  see  Figure  3- 

As  we  have  pointed  out  on  several  occasions,  the  Navier- 
Stoke  equations  may  be  extracted  from  the  thirteen  moments 
equations.   This  may  be  done  by  direct  inversion  of  the 
system  corresponding  to  (3.19)-   In  that  case,  the  point  of 
view  is  that  one  obtains  a  system  of  differential  equations 
which  gives  the  hydrodynamical  branch  accurately  to   0(k  ). 
More  generally  (since  it  can  be  applied  to  non-linear  systems) 
one  can  obtain  it  by  the  interpolation  scheme  given  by  Grad. 
(See  reference  3-  P-  110).   In  any  case  the  system  of 
equations  is  given  by, 


St  +  3x 


u 


0 
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(5-17)    -+Bf  (P+T)  +5 


4 afu 

A02   Sx2 


=  0 
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bT        2      du 
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The  dispersion  relation  of  the  Navier-Stokes  systems  is  of 
course  given  by  (3. 19).   It  contains  only  the  hydrodynamical 
branch,  which  it  gives  faithfully  for  small  wave  number.   The 
dispersion  relation  sketcked  for  all  wave  number  is  given  in 
Figure  4. 

The  behavior  near  the  origin  is  given  to   0(k  )   by 
(3.23)  and  (3.24).   A  simple  calculation  shows  that  for 
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large  k  ,   the  three  roots  are, 

(5-18.)    a-  =    D\\^ +  0(l/k2) 

DA02      11 


(5.19)  °  =  ^T"  +0(1) 

^02 

(5.20)  o  =  ^L  +   0(1) 

;A11 


So,  as  is  sketched  In  Figure  4,  the  two  propagating  modes 
return  to  the  real  axis  in  a  double  point  and  proceed 
independently  along  the  real  axis.   This  behavior  could 
also  have  been  noticed  earlier  from  equations  (5.17).   There, 

if  k  >  00   (or  equivalently  -^  )  <*>)   the  single  space 

derivatives  in  the  last  two  equations  are  ignorable  with 
respect  to  the  dissipative  derivatives.   And  the  phenomena 
is  essentially  governed  by  two  diffusion  equations. 

This  state  of  affairs  is  somewhat  ironic.   The  Navier- 
Stokes  equations  have  been  extracted  from  the  thirteen 
moments  equations.   On  comparison  of  Figure  3  and  4,  the 
Navier  Stokes  equations  give  a  bad  approximation  to  the 
latter  equations  for  high  frequency  phenomena.   They,  in 
contrast  to  the  moments  equation,  give  unbounded  absorption 
with  increasing  wavenumber.   On  the  other  hand  the  moments 
equations  give  a  bad  approximation  to  the  Boltzmann  e- 
quation  for  large   k  .   The  irony  of  course  is  that  the 
Navier-Stokes  equations  qualitatively  at  least  approximate 
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the  Boltzmann  equation,  for  large  k. 

The  remarks  made  thus  far  indicate  how  the  thirteen 
moments  might  be  improved  so  as  to  give  a  more  faithful 
approximation  to  the  exact  equations.   First  it  is  obvious 
that  the  relative  placing  of  the  eigenvalues  must  be  given 
consideration.   Prom  relaxation  considerations  it  is  clear 
that  in  a  homogeneous  problem,  the  moment   bpo   persists 
for  a  longer  time  than  does  the  stress   b    (since  it  is 
represented  by  a  smaller  eigenvalue).   Looking  at  Table  II 
the  eigenvalue  ordering,  it  is  tempting  to  consider  the 
truncated  system  i,j  <  2.   This  for  small   k  would  give 
the  first  two  clusters  of  branches  in  Figure  1  .   The 
immediate  shortcoming  to  this  system  is  that  it  excludes 
the  stress.   This  would  have  the  consequence  of  giving  the 
hydrodynamic  branches  incorrectly  to   0(k  ).   Alternately 
the  Navier-Stokes  could  not  be  extracted  from  this  system. 
A  more  plausible  system  which  rectifies  this  difficulty, 
is  given  by  the  first  four  equations  of  (5.1)  and  the 
following  two, 

(5-21)    ^  -  AniS,  +  ^  +  §  |T+!_0 
v     '     dt     11  1     ox     2   ox    ox 


(5.22)     |^Q  -  A2QQ+i  ^Sl=  0 


where 
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(5-23)        Q  =   )/lO/5  b 


20 


(2tt) 


1        -/e-^1^ 


372- 


5I2+  |5] 


The  system  given  by  (5. 21),  (5.22)  and  the  first  four 
equations  of  (5.1),  still  shares  with  all  truncated  systems 
the  shortcoming  of  having  bounded  absorption.   One  way  to 
overcome  this  is  to  use  the  interpolation  scheme  suggested 
in  the  extraction  of  the  Navier-Stokes  from  the  thirteen 
moments  equations.   In  place  of  the  fourth  equation  of  (5.1) 
we  take  instead 


(5-24)    p 


4   du 


as, 


11  '    3AQ2  "3x   15^02  "5x~ 


Combining,  the  entire  system  reads, 


s 


(5.25)  < 


!t.+&-° 

lr  +  a|(p+T)  - 

4 

^2 
d   u 

>    2 

8 

a2sx 

15?w 

>  2 
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ST    2   du    2   aSl 
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-  AnnS 


^2 
d  u 


d2S 


11  !  '  ^02  ax2   "  ^^2     Sx2 


1     5    ^t     ao 

"  +  I  ^x"  +  ^x" 


=  0 


-  A._Q  +  ^ 


4     a 


*20 


3"  "Sx"  1 
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This  system  is  easily  generizable  to  a  non-linear  system. 

We  will  not  study  this  nor  will  we  study  (5.25),  in  this  report, 
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We  are  interested  in  the  more  general  approximations  to  the 
Boltzmann  equations  which  will  be  made  shortly. 

Implicit  in  each  of  the  truncation  and  associated  interpo- 
lation methods  is  a  particular  form  of  the  distribution  function. 
Equivalent  to  the  method  of  straightforward  truncation  is  the 

choice  of  a  distribution  function  of  the  form, 
N 

(5.26)     g  =  3~  b  if 
r=o 

where  some  particular  ordering  is  taken.   This  is  substituted 

into  the  linearized  Boltzmann  equation,  and  moments  are 

24 
taken.   One  then  finds  that  the  first   (N+l)     equations 

form  a  determined  system  for  the  coefficients 

b^(r  =  0,1,  ...;N)   which  is  identical  to  the  truncated  system. 
After  solving  this  system  of  equations  for  the  b's,  they  are 
substituted  into  (5.26)  above,  and  the  corresponding  distri- 
bution function  is  known.   This  is  in  fact  the  general  pro- 

p 
cedure  used  by  H.  Grad  .   In  order  to  obtain  the  interpolation 

scheme,  one  takes  the  interpolation  equations  and  substitute 

them  into  equation  (5.26).   We  illustrate  this  procedure  for 

both  the  Navier-Stokes  equations  and  the  thirteen  moments 

equations . 


24.    Actually  more  than   (N+l)   equations  may  be  generated, 
The  embarrassing  situation  of  having  more  equations 
than  unknown  is  obtained.   This  can  be  avoided  by 
ignoring  equations  without  time  derivatives. 
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In  order  to  obtain  the  thirteen  moment  equations,  one 

2 

takes,  with  H.  Grad   ,  for  the  distribution  function, 


(5-27)    g  =  P  +  u1^1  +  2^   (-5-  -  1) 


*Pll(ej  -V  +  8^^- 1) 


This  is  substituted  into  the  linearized  equation,  moments 
are  taken  and  the  system  (5«l)  is  obtained.   On  solving  the 
latter  system  one  substitutes  the  solution  into  (5.26)  and 
the  distribution  in  this  approximation  is  known..  As  pointed 
out  earlier  in  this  section,  in  the  Navier  Stokes  interpolation 
one  uses, 


<^8)   Pn^i 


(5-29)    \  =  it^~%L 


These  latter  two  equations  are  substituted  into  (5.26),  to 
get 


/\  /\       f* 


'   %,*?    .  -m  ^  1  .du 


(5O0)     g  =  p  +  »£,  +  2.   Tp-  -  1)  +  -i_  2± 

1  1  2    3  A02  ^ 

^1     3   An  Sx  ^  2    2^ 

This  is  substituted  Into  the  linearized  Boltzmann  equation, 
moments  are  taken  and  the  Navier  Stokes  equations  result.   The 
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latter  procedure  is  analogous  to  the  Chapman-Enskog  procedure. 

The  remarks  made  about  the  truncated  systems  now  apply 
to  the  corresponding  distribution  functions.   Essentially 
that  they  can  only  be  pressumed  to  be  asymptotically  valid. 
One  further  remark  in  support  of  this  is  the  following. 
In  demanding  that  a  distribution  be  of  the  form  (5.27)  or 
(5.30)  we  eliminate  consideration  of  free  flow  or  near  free 
flow,   for  the  insertion  of  such  a  distribution  function  in 
a  free  flow  operator  immediately  leads  to  the  generation  of 
the  higher  moments  which  have  been  excluded  from  consideration, 

We  now  investigate  another  type  of  approximation  to  the 
Boltzmann  equation  which  is:  (l)   asymptotically  correct  ', 
(2)   does  take  into  consideration  the  higher  moments  ',       (3) 
and  gives  a  transition  to  free  flow. 

We  begin  by  considering  a  relaxation  theory,  i.  e.  the 
space  homogeneous  case.   Setting  the  space  derivatives  of 
(2.4l)  to  zero  we  get 


a 


(501)    (^fr  -  a  >nv  =  0 


"ST" 


M-V 


Since  all  but  the  first  three  A's  (which  are  zero)  are  nega- 
tive the  moments   b    vanish  exponentially.   Since  |A   I— \   00 

U.V  a  ■  U.V '     ' 

the  higher  moments  vanish  more  rapidly.   It  seems  plausible 

therefore  if,in  the  space  inhomogeneous  case  we  choose  all 

A  .  .   after  some  point  to  be  a  single  constant,  that  we  will 
u. '  v  '  ^  to 
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overestimate  the  effect  of  the  higher  moment.   The  relaxation 
theory  would  then  be  correct  for  b   's   which  occur  before 

U.V 

b  ,  .  .   For  those  which  come  after,  it  is  qualitatively 
correct,  giving  exponential  decay,  but  with  the  incorrect 
decay  time. 

We  carry  out  this  operation  in  connection  with  equation 
(2.36).   Instead  of  the  latter  equation  we  write, 


(5.32)     (-A-  +  ?n  3^7)  y_  b  of    0   =  2_   b  of    0^ 

K^    ^      '  Mt1      1  OX  '  '    *—n        M  M    ««- MM  M 

r,  i!  0,0 


r,i>p,q 


The  remarks  we  are  now  making  are  for  the  polynomial  order, 
and  under  this  manipulation,  Table  III  changes  only  in  the 
diagonal.   After  some  points  on  the  diagonal  the   A   »  '  s 
are  replaced  by   A  ,   which  is  still  to  be  specified.   We 
take   A  to  be  the  next  eigenvalue  so  as  to  accurately  give 
the  equation  for  one  more  moment.   The  advantage  to  equation 
(5.32)  lies  in  a  form  into  which  it  may  be  transformed. 
Adding  "zero"  to  the  right  hand  side  of  (5.32),  we  have, 

pq  P,q 

a       .  />    a- 

Jr^ri 

00  <XI 


(5-33)  (^It  +  ^  3§)g  -  }_  \ifvSi\SL   ~   *  2-  \if 


c.-.i 


p,q 


0,0  00 
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Equation  (5-33)  is  precisely  the  equation  developed  by  Gross 

6 
and  Jackson  as  a  generalization  of  the  equation  found  in  the 

4 
paper  by  Bhatnager,  Gross,  and  Krook  .   As  an  example  of 

(5- 33) ,    we  take   A  =  A„   and   (p,q)  =  (1,0)  ,   then  we  get 

~2 
(5-34)   (^  +  ^  J  -  AQ2)g  =  -  AQ2[p  +  u1i1  +  |  T(^  -  1)] 


This  is,  the  linearized  form  of  the  single  relaxation  model 
found  in  references  4  and  6, 


(5-35)     (f§  +  «! ,^)f  =  -A02(V  f ) 


where   f   is  the  local  maxwellian, 
o  ' 


-U-u(x)]2 


/   ,  \      2RT(x 
(5-36)  f  =    p{X}tK/2   e 
0    [2ttRT(x)]3/2 


The  linearization  is  that  given  by  (2.4)  and   -Afi?  is  the 
"collision  frequency".   Another  model  worth  noting  is  gotten 
by  taking   A  =  A    and  (p,q)=  (l,l),   then 

(5-37)    (^|r  +  *'  J  -  AQ3)g  =  -A^tp+f^  |  Tt-L  -  1) 

A  /n    »2        A-,-,^/\£2 
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Appropriate  names  for  (5-35)  and  (5-37)  might  be  the  Euler 
and  thirteen  moments  models  respectively.   The  first  model 
(5- 3^) ?  for  the  problem  of  sound  propagation,  was  studied 
in  reference  (4).   A  much  more  detailed  analysis  of  sound 
propagation  using  this  model  has  been  performed  by  the 
author  and  J.  Thurber  ■*    .   In  that  paper  (5-37)  is  also 
applied  to  sound  propagation. 

We  recall  that  these  models  (5-35)  and  (5-37)>  were 
obtained  from  Table  III  by  replacing  all   A's   past  a  certain 
point  by  a  single  constant.   Let  us  consider  the  hydrodynamic 
branch  of  these  models.   First  we  consider  equation  (5-35)- 
Referring  to  Table  III,  we  see  that   i,j  <  1   remains  unchanged 
and  the  three  branches  in  the  hydrodynamic  cluster  is  given 
correctly  to  0(k).   To  the  next  order  however  we  see  that 

and 

the  expansions  for   o  are  somewhat  modified.  In  fact  all 
higher  orders  are  modified  in  this  way.   However   the 
complexion  of  the  matrix  of  Table  III  has  not  been  changed, 
and  the  series  given  by  the  single  relaxation  model  is  of 
the  correct  form  but  with  incorrect  constants.   A  truncation 
theory  for  this  case  may  be  developed  as  in  the  case  of  the 

25.  L.  Sirovich  and  J.  Thurber,  "Sound  Propagation  according 
to  Kinetic  Models".  NYO-9757,  Institute  of  Mathematical 
Sciences,  New  York  University  (August,  1961). 

26.  There  is  some  arbitrariness  at  this  point.   We  have 
chosen   A  in  (5.32)  so  that  the  propagating  branches 
are  correct  to   0(k2).   In  this  case  the  diffusing  mode 
is  given  inaccurately.   Equivalently  we  could  give  the 
diffusion  and  inaccurately  give  the  propagation.   The 
pertinent  point  is  that  we  have  just  one  parameter  left 
open  to  us.   We  can  correctly  specify  the  viscosity  of 
the  heat  conductivity,  but  not  both  in  general. 
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A,,   is  replaced  by   An?   and  the  expansions  for   A0„ 


exact  equation.   The  lowest  truncation  1-,  j  <  1   gives,  of 
course,  the  Euler  dispersion  relation  of  Figure  1  .   This  is 
also  indicated  in  Figure  5,  where  the  effect  of  successively 
higher  truncations  on  the  hydrodynamic  branches  are  sketched 
by  the  doited  lines  moving  to  the  left.   The  diffusion  mode 
moves  to  the  left  along  the  real  axis,  as  is  shown.   The 

exact  theory  of  the  single  relaxation  model  is  given  by  the 

27 
heavy  lines  terminating  in  Ra  =  A„   ' .   From  the  remaks 

made  earlier,  we  know  that  no  truncated  system  can  have  a 
greater  absorption  than  that  given  by   Aop.   For  this  reason 
that  line   Rea  =  ~h  plays  the  role  it  does  in  Figure  5  . 

Nothing  has  been  said  of  the  infinity  of  branches  which  have 
their  origin  at   a  =  Anp  .   For  any  finite  truncation  these 
extraordinary  branches  appear  in  the  strip   A0?  <  Ra  <_  0  . 
As  the  order  of  the  truncation  increases  these  are  shoved 
over  to  line,   Ra  =  A   ,   and  in  fact  the  latter  line  be- 
longs to  the  spectrum. 

It  is  also  revealing  to  consider  the  systems  of  partial 
differential  equations  gotten  by  truncation  and  interpolation, 
as  was  done  for  the  exact  Boltzmann  equation.   As  was  indicated 
the  Euler  equations  drop  out  as  the  first  truncation.   The 
thirteen  moments  come  out  as  the  second  truncation,  however 
with  only  the  one  constant.   For  the  first  interpolation 
the  Navier-Stokes  result,  however  again  with  just  one  constant. 


27.   This  line,   Rea  =  ~h        ,      is  also  included  in  the  spectrum, 
see  reference  25- 
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One  may  choose  it  to  agree  either  with  the  viscosity  or  the 

28 
heat  conductivity.   This  has  been  already  remarked  upon 

One  gets  the  higher  moment  equation  and  the  higher  inter- 
polation correctly,  but  with  the  remark  about  one  constant 
always  applying. 

We  can  comment  in  the  same  way  about  equation  (5- 37)- 
Again  referring  to  Table  III,  all  A's   after   A    are  set 
equal  to   An^  .   Recalling  the  calculations  made  in  connection 
with  the  exact  theory  we  see  that  the  Euler,  Thirteen  Moments, 
Navier-Stokes,  and  Burnett  equations   of  equation  (5-37)  are 
all  identical  with  the  similar  forms  gotten  from  the  exact 
Boltzmann  equation.   The  comparable  statement  can  be  made  in 
reference  to  the  dispersion  relation.   The  hydrodynamical 
branches  are  given  correctly  to   0(k  ).   A  sketch  of  the 
dispersion  relation  is  given  in  Figure  6  .   One  can  identify 
the  first  set  of  dotted  lined  with  the  thirteen  moments 
dispersion  law  given  in  Figure  3  .   Successively  higher 
approximations  are  given  by  the  dotted  lines  which  go  to  the 
left.   The  heavy  dark  line  is  the  exact  spectrum  sketched 
from  the  results  in  reference  (25).   Again  the  peculiar 
effect  of  the  circled  region  in  Figure  3,  is  exhibited.   The 
effect  of  taking  more  moments  is  to  move  this  double  point 
to  the  left.   The  same  sort  of  remarks  can  be  made  now,  as 
in  the  previous  case,  and  again  the  reader  is  referred  to 
reference  (25)  for  a  detailed  investigation. 


28.   L.  Sirovich  "On  the  Kinetic  Theory  of  Steady  Gas  Flows" 
NYO-9086,  Institute  of  Mathematical  Sciences,  New  York 
University  (i960). 
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Once  again  the  shortcomings  of  choosing  a  polynomial 
ordering  occur  in  connection  with  the  branches  starting  at 
A-,-,   and  A„   .   We  now  offer  a  method  alternate  to  the  one 
suggested  by  Gross  and  Jackson  . 

The  discussion  parallels  that  used  in  modifying  the 
thirteen  moments  equations.   From  Figure  1  we  see  that 
equation  (5-37)  does  not  represent  the  exact  dispersion  law 
in  the  neighborhood  of   A-,,  ,  Afi   .   In  order  to  rectify  this 
we  choose  the  ordering  implied  by  the  eigenvalues  according 
to  Table  II.   Instead  of  equation  (5.37)  we  want  to  include 
the  moment  Q  as  defined  by  equation  (5.23).   Mimicking  the 
arguments  used  in  developing  the  equations  given  by  (5. 21), 
(5-22),  we  take  our  model  equation, 

(5-38)       (3!  +  ^  a|  -  A02)g  -  -a02[p~  +  l±  u±  + 1  J(£  -  1) 

(l  -  ^ii)si  (4  -  1)  +  (1  -  ^2)Q(ip  -  5I2+  ^f) 
A02        D  02 

Referring  to  Table  II  we  can  make  the  following  statements 
about  equation  (5.38).   On  setting  the  A's   from   Anp   onward 
equal  to  ~hQ        we  only  change  the  diagonal  for   i,j  >  3  • 
That  is,  the  relaxation  times  of  eight  moments  (ij  <  3)   is 
given  faithfully.   This  is  to  be  compared  with  (5.37)  which 

gives  six  relaxation  times  correctly.   Furthermore  the 

3 

hydrodynamic  branches  are  given  correctly  to   0(k  )  with  the 

accuracy  as  given  by  (5.37).   On  the  other  hand  the  second 
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branch  is  given  to   0(k)  by  (5.38),  but  not  at  all  by  (5«37)> 
as  already  mentioned.   In  connection  with  the  truncated  and 
interpolated  equations,  we  can  obtain  the  Euler,  Navier- 
Stokes,  Burnett  and  thirteem  moments  systems  from  (5-38),  all 
which  are  also  obtainable  from  (5. 37  )•   However  (5-38)  also 
furnishes  the  system  given  in  connection  with  (5.21),  (5-22) 
and  also  the  associated  interpolated  system  given  by  (5-25)- 
Although  the  dispersion  relation  associated  with  (5.38)  has 
not  been  studied  yet,  it  is  reasonable  to  conclude  that  it 
will  have  a  dispersion  relation  of  the  form  sketched  in  Figure 

7- 

It  is  interesting  to  note  that  both  (5-37)  and  (5-38) 
have  the  same  number  of  moments  on  the  right  hand  side. 
However  in  order  to  rectify  the  same  difficulty  in  connection 
with  the  thirteen  moments  equations  it  was  necessary  to 

introduce  another  moment  Q.   This  simplification  is  due  to  the 

6 
device  of  Gross  and  Jackson  ,  used  for  equation  (5. 33)-   This 

allows  us  to  include  accurately  at  least  one  more  moment  than 
is  given  explicitly  on  the  right  hand  side  of  the  model 
equation. 

It  should  not  be  construed,  that  for  purposes  of  trun- 
cation, interpolation  or  constructing  models,  that  the  eigen- 
value ordering  is  any  more  suitable  than  the  polynomial 
ordering  which  is  generally  used.   For  instance,  in  almost 
any  conceivable  problem  in  fluid  mechanics  the  hydrodynamic 
mode  is  all  one  is  interested  in.   As  we  saw  in  Section  3, 
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the  polynomial  ordering  is  the  most  natural  for  determining 
the  hydrodynamic  branches.   An  entirely  different  ordering 
(neither  eigenvalue  or  polynomial)  would  be  used  if  one 
wanted  to  examine  another  branch.   Examples  of  this  are 
given  by  Tables  IV  and  V.   In  any  problem  for  which  either  a 
model  equation  or  a  truncated  system  is  to  be  used,  one  must 
decide  how  much  of  each  branch  to  describe  exactly.   Going 
back  to  our  discussion  of  the  thirteen  moments  we  can  make 
this  more  concrete.   There  we  saw  that  the  thirteen  moments 
truncation  is  the  one  necessary  for  the  exact  description  of 
the  hydrodynamic- branches  to   0(k  ).   In  the  process  however 
the  circled  area  in  Figure  3   appears.   This  seriously  distorts 
the  same  area  in  the  exact  picture  given  by  Figure  1  .   There 
were  then  two  immediate  courses  open  to  us  if  we  wished  to 
rectify  this  situation.   We  could  drop  the  stress  and  stress 
equation  and  replace  it  by  Q  and  the  associated  equation. 
This  introduced  no  more  complexity  in  as  much  the  number  of 
equations  remained   equal  to  five.   However  it  reduced  the 
accuracy  of  the  hydrodynamic  branches  to   0(k  ).   If  by 
prearranged  agreement,  no  more  than  five  equations  are  allowed, 
and  accuracy  in  the  hydrodynamic  region  is  required  then  only 
the  thirteen  moments  system  must  be  taken.   We  must  then  be 
content  with  a  poor  description  of  the  circled  region  of 
Figure  j5  •   It  seems  a  moot  point  as  to  whether  a  poor  de- 
scription is  better  than  no  description  at  all.   In  the 
latter  case  we  might  as  well  discard  the  thirteen  moments 
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system  in  favor  of  an  interpolation  such  as  Navier-Stokes  or 
Burnett.   If  on  the  other  hand  we  allow  ourselves  six 
equations,  we  can  use  the  system  suggested  earlier  and 
connected  with  equation  (5.22)  or  its  associated  interpolation 
(5«25).   Each  is  gotten  from  the  eigenvalue  ordering  of  Table 
II  (each  of  these  gives  the  second  cluster  to  0(k).   The 
hydrodynamic  roots  in  the  next  truncation  are  gotten  to   C^k-5) 
which  is  the  system  of  eight  equations  comparable  to  the 
8X8   truncation  of  Table  III. 

In  general  there  is  no  systematic  way  in  which  to  order 
the  eigenf unctions.   This  can  be  more  clearly  pointed  out  If 
we  allow  ourselves  the  luxury  of  a  large  number  of  moments  in 
our  description.   Suppose  for  example  we  choose  to  use  the 
polynomial  ordering.   We  would  have  then  the  hydrodynamical 
branch  given  quite  accurately  but  we  would  have  also  the 
incorrect  evolution  of  a  great  number  of  moments.   Many  of  the 
moments  which  occur  between  those  involved  in  the  hydrodynamical 
calculation  will  not  even  appear.   On  the  other  hand  strict  use 
of  the  eigenvalue  ordering  leads  to  a  correct  dispersion  for 
samll  k  ,   but  does  not  necessarily  give  any  branch  to  a  great 
accuracy.   A  compromise  between  the  two  methods  might  be 
advantageous.   One  choses  which  branches  to  describe  with  some 
degree  of  accuracy  and  then  collects  the  moments  accordingly. 
This  of  course  corresponds  to  the  neither  the  eigenvalue  nor 
polynomial  ordering. 

We  close  by  considering  the  model  equations  from  a  formal 
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point  of  view.   The  following  remarks  are  independent  of  the 
ordering  chosen.   This  being  the  case  we  represent  a  model 
equation  "by, 

N 

(5-59)   (J  +  ^  +  v)g  =  $I  VrAi 

n=0 

The  identification  with  the  previous  equations  are  easily  made, 
vn  being  related  to  (l  -  -rip),   and  i/        being  the  eigen- 
functions   ^  „    .      If  we  regard  the  coefficients   b   as  known 
functions  of   (x,t),  equation  (5-39)  can  be  easily  integrated. 
The  result  is 

t 

(5-40)    g(x,t,|)  =  e~vt  J  evs  >   b  (x+Js-t),s)lv  ds 

o      n=0         ^   J 


+  g(x  -  4^,1,0) 


-vt 


the  last  term  being  the  initial  data.   As  the  integrated 
equation  indicates  the  integration  is  along  particle  paths, 
We  can  now  operate  on  (5.40)  with 


(5-^1)    /  ^nf°  d|    for    n  =  0,1,  •••,  N 

to  get  a  system  of  linear  equations  in  the   N  +  1   variables 


29 


bn(x,t),  n  =  0,1,  .  .  .  N 

29.    This  may  be  reduced  to  a  system  of  Abel  integral 
equations. 
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On  solution  of  the  integral  equations,  the  quantities  are 
substituted  into  (5.40)  and  the  perturbed  distribution  function 
is  obtained.-   We  illustrate  the  procedure  for  N  =  2,   the 
single  relaxation  model  (5.;4).   Letting   v  =  -A^,  the 
solution  given  by  (5-40)  is 


vt  f  vs 


(5-^2)   g(x,t,0  =  e     Je   [p(x+l1(s-t|  ,s)+^1u(x  +  ^1|s-tj  ,s) 


+  T(x+l1|s-t],s)(-^2  -  5)]ds 


+  g(x-£.jt,€  0)e 


■vt 


Then  operating  on  (5-42)  with  (5-4l)  for  n  =  0,1,2,  we  get 
the  system 


(5-43) 


p(x,t) 
u(x,t) 
T(x,t) 


■vt 


J 


■  vt 


+  e 


vs 


^  2 


sl 
(|-  3) 


g(x-e,t,^10)f°de 


'1 


-  3)  ^(- 


.2 


(^  "  ^ 
P2 


3)2  (^-3)2 


p(s) 

u(s) 

T(s) 


ds 


where  the  argument  of  the  elements  in  the  last  vector  is 
(x+£u{s-tj  ,s  ).   For  instance   p(s)  =  p  (x+^-J  s-tj  ,  s  ) . 
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One  solution  of  this  system,  the  flow   (p,u,T)  is  substituted 
into  (5- ^2)  giving  the  perturbed  distribution  function.   We 
see  that  the  distribution  function  is  given  in  terms  of  the 
initial  distribution  and  the  evolution  of  the  flow   (p,u,T). 
In  general  any  model  gives  the  distribution  function  in 
terms  of  the  initial  distribution  and  the  evolution  of  a 
finite  number  of  moments.   In  a  certain  sense  this  is  a 
generalization  of  the  moments  method  of  Grad.   There,  as  we 
saw  earlier, ( 5. 27) ,    the  distribution  is  given  in  terms  of  a 
finite  number  of  moments  and  one  solves  differential  equations 
for  the  moments.   However  full  use  of  the  initial  data  is  not 
made  there,  and  as  our  considerations  indicated,  the 
distribution  function  is  only  asymptotically  valid.   On  the 
other  hand  the  model  solution  describres  the  flow  as  we 
approach  free  flow.   In  fact  on  setting  the   v's   to  zero 
in  (5-^0)  we  get  the  free  flow  solution. 
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APPENDIX  I :   Two  Lemmas  on  Determinants 

Notation 
Let 

(1.1)  D°  =  D(a   ) 

be  an  n  x  n  determinant  with  elements   a    which  are 
dif ferentiable  functions.   We  write 

da 

(1.2)  -^  =  a'    . 

v    '  dx       |iv 


By 


°(<v.) 


we  shall  mean  the  determinant  with  the  \xv      element  replaced 
by  its  derivative.   Further  we  write 

n 

(1-3)     D1  =  X  D(a'  ) 

u.v 

where   the    summation   is    over  all   elements.      In  general, 

n 
(I  A)  Dk   =     X  D(Vv/     •••'    Vv  J 

where  the  summation  is  over  all  values  of  each  pair  of  sub- 
scripts.  If  an  element  appears  m   (m  <L  k)   times,  it  denotes 
the   m   derivative.   With  this  notation  our  fitst  lemma  is; 
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Lemma   1 . 

(1-5) 

dm 
dxm 

D  =  (D1-  n(n-l))mD° 


The  right  hand  side  of  (l-5)  represents  the  formal  binomial 
expansion,  with  D   given  by  (l.4). 

To  prove  this  we  start  with  the  well-known  expansion 
for  the  derivative  of  a  determinant 


n 
hD   =  5  D(a^,  a^,  ...,  a^n; 


but 

n 


i 

.a. 


D(a^,  a^,  ...,  a^)  +  (n-l)D°  =  |~  D(«  _. 

and  on  summing  this  over  all  rows  we  get 

(1.6)      ^  D  =  D(l)  =  D1  -  n(n-l)D°  =  [D1-  (n-l)n]D°  . 

Differentiating  (1.6)  and  applying  it  to  itself  we  see  that 
(1.5)  is  true  for  m=2.   Using  the  inductive  hypothesis,  we 
have 


dnD   n(n)    dn  1      j.1        r_,l    /   ,nn-l  nl 
— n  =  D    =  — ^ZT   D  =  [D  -  n(n-l)]     D 
dx  dx 


=  [D1-  n(n-l)]n-1  (D1-  n(n-l))DC 


=    [D1-  n(n-l)]n  D° 
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and  the  lemma  is  proven. 

In  our  work  we  are  interested  in  evaluating  the  determi- 
nant  D(a,k)   and  its  derivatives  at   k=0.   Further,  according 
to  equation  (3-13)  we  get 


(1.7)      0  =  D°  =  D1  =  D2  =  ...  =  Dn 


Our  second  lemma  gives  a  sufficient  condition  for  a 
determinant  to  yield  an  equation  with  real  coefficients.   We 
consider  determinants  whose  elements  are  either  real  or 
Imaginary.   A  zero  element  is  considered  as  having  both 
parities . 

Lemma  2 .   A  determinant  for  which  the  parities  of  a  row  may 
be  gotten  by  multiplying  the  preceeding  row  by  I ;  or  any 
determinant  which  may  be  put  into  this  form  (say  by  extracting 
i's  from  rows)  has  one  parity.   That  is,  the  expansion  of  the 
determinant  results  in  a  form  all  of  whose  coefficients  have 
the  same  parity. 

The  lemma  is  clearly  true  for  a  two  by  two  determinant. 
Expanding  an  n  x  n  determinant  (which  staisfies  the  hypo- 
thesis) along  the  first  column  we  get   n   (n-l)   (n-l) 
determinants.   The  first  determinant  of  the  expansion  satisfies 
the  hypothesis,  the  second  determinant  can  be  given  the  same 
form  by  interchanging  the  first  and  second  rows.   In  fact  each 
determinant  of  the  expansion  can  be  put  into  the  form  of  the 
first  determinant  by  interchanging  rows  and  extracting  an  I 
when  necessary.   Imposing  the  inductive  hypothesis  we  prove  the 
lemma. 
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Figure  2.    THE   NON-DISSIPATI VE  OR 
EULER   DISPERSION   LAW 


Figure  3.    THE   THIRTEEN   MOM 
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